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Abstract 

The  relationship  between  the  observed  increase  in  the 
flux  of  low  energy  (<500  ev)  electrons,  at  an  altitude  of 
800  km  and  the  corresponding  amplitude  and  frequency  scintil¬ 
lations  of  radio  waves  passing  through  ionospheric  irregular¬ 
ities  at  350  km  was  studied.  This  data  was  measured  by  the 
electron  spectrometer  (J-sensor),  the  magnetometer,  and  the 
radio  beacon  experiment  (137  MHz  and  413  MHz)  of  the  HILAT 
satellite.  The  data  analyzed  here  was  taken  from  January- 
March  1984.  During  that  time  interval,  198  enhancements  of 
electron  flux  were  recorded  in  the  night  sectors  (magnetic 
local  time),  occurring  most  frequently  at  0300  and  2300  MLT . 

For  73  of  these  enhancements,  magnetometer  data  was  exam¬ 
ined  for  corresponding  geomagnetic  disturbances  (dBy/dt),  but 
only  nine  had  such  associated  disturbances.  The  average 

increase  in  electron  number  flux  for  the  73  enhancements  was 

8  2 

found  to  be  9.7  x  10  electrons/cm  -sec-ster  with  a  standard 

8  2 

deviation  of  9.6  x  10  electrons/cm  -sec-ster.  Enhancements 

9  2 

with  number  fluxes  greater  than  2  x  10  electrons/cm  -sec-ster 
were  the  same  ones  that  had  correspondingly  large  dBy/dt 
values . 

When  the  values  of  the  phase  scintillations  for  137  MHz 
(known  as  PRMS)  were  plotted  against  those  for  413  MHz,  a 
linear  curve  that  agreed  with  the  theory  put  forth  by  Briggs 
and  Parkin  resulted.  The  same  was  done  with  the  amplitude 


scintillation  data;  the  distribution  was  not  as  linear,  but 
many  points  fell  along  linear  curves  for  limiting  large  and 
small  values  of  the  traverse  radii  of  irregularities  as  de¬ 
scribed  by  Briggs  and  Parkin. 

For  the  42  sets  of  values  for  PRMS  and  S4  used  in  the  dis¬ 
criminant  analysis,  the  increase  in  integrated  number  flux 
(PEAK)  and  the  average  electron  energy  (EAVG)  were  significant 
indicators  in  determining  S4  at  137  MHz.  The  size  of  the 
enhancement,  as  indicated  by  the  time  necessary  to  traverse  it 
(TIME)  and  the  time  difference  between  enhancement  encounter 
and  scintillation  occurrence  (DT)  were  good  indicators  for 
determining  the  rise  of  PRMS  at  both  frequencies  and  the  rise 
of  PRMS  at  413  MHz  by  itself. 

TIME  and  DT  were  also  indicative  of  scintillation  increases 


(S4  and  PRMS)  at  413  MHz.  Results  for  finding  indicators  for 
S4  at  413  MHz,  and  PRMS  at  137  MHz  were  not  statistically  sig- 
nif icant . 


STATISTICAL  ANALYSIS  OF  LOW  ENERGY  ELECTRONS 


AND  THEIR  RELATIONSHIP  TO  RADIOWAVE 
SCINTILLATIONS 

I .  Introduct ion 

Background 

In  order  to  deter  a  nuclear  attack  against  the  United 
States  it  is  essential  to  have  a  credible  capacity  to  retali¬ 
ate  against  such  an  attack  from  any  potential  enemy.  One 
necessary  ingredient  for  maintaining  such  a  force  are  command, 
control,  and  communications  (C3)  systems  that  can  reliably 
function  during  such  a  conflict  and  its  aftermath  (Fremouw, 
1983:3) . 

A  major  source  of  degradation  of  our  C3  systems  would  be 
radio  scintillations  associated  with  high  altitude  ionospheric 
plasma  irregularities  created  by  nuclear  bursts.  Since  the 
Defense  Nuclear  Agency  (DNA)  is  charged  with  investigating 
such  adverse  nuclear  effects,  the  development  of  these  plasma 
irregularities  are  of  major  interest  (Fremouw,  1983:5). 

Since  atmospheric  nuclear  explosions  are  forbidden  by  the 
Nuclear  Test  Ban  Treaty,  it  is  indeed  fortuitous  that  such 
irregularities  occur  naturally  under  solar  produced  disturbed 
geophysical  conditions.  Although  much  active  experimentation 
and  observation  has  been  done  on  these  irregularities,  much 
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remains  to  be  learned  about  how  these  plasma  structures  develop, 
evolve,  and  decay.  In  fact,  these  processes  are  especially 
more  complicated  at  the  higher  latitudes  than  at  the  equatorial 
or  middle  latitudes.  This  is  mainly  due  to  coupling  that 
occurs  between  the  ionosphere  and  the  magnetosphere.  For  in¬ 
stance,  it  is  believed  that  precipitating  electrons  lead  to 
the  establishment  of  the  scintillation-producing  irregulari¬ 
ties,  yet  the  magnetospheric  processes  responsible  for  this 
precipitation  are  poorly  understood  (Rino,  1983:1179). 

Since  the  Air  Force  Geophysics  Laboratory  (AFGL)  is  tasked 
with  the  study  of  geophysical  phenomena  and  how  such  phenomena 
affect  Air  Force  systems,  it  would  seem  then  that  a  joint 
investigation  of  radio  scintillations  produced  by  enhanced 
plasma  structures  would  be  a  well-suited  endeavor  for  both 
AFGL  and  DNA. 

The  DNA  Wideband  Satellite  Experiment  launched  in  1976 
transmitted  several  phase-coherent  radio  signals  which  provided 
useful  three-dimensional  information  on  the  configuration  of 
these  irregularities.  However,  in  order  to  perform  the  pre¬ 
viously  mentioned  study,  an  orbiting  radio  propagation  experi¬ 
ment,  along  with  other  on-board  instrumentation,  is  needed  to 
explore  the  plasma  and  free  energy  sources  that  are  associated 
with  ionospheric  irregularities  (Fremouw,  1983:4). 

The  response  to  this  requirement  is  the  DNA- sponsored 
satellite  known  as  HILAT  (High  Latitude).  HILAT’s  compliment 
of  experiments  consists  of  a  multifrequency  radio  beacon,  an 
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energetic  particle  detector,  a  magnetometer,  a  plasma  probe, 
a  driftmeter,  and  an  ultraviolet  imager.  Hopefully  this  will 
provide  data  for  a  more  complete  description  and  understand¬ 
ing  of  high  latitude  plasma  structures  than  before  possible. 

The  characterization  of  the  complex  radio  signal  (ampli¬ 
tude  and  phase)  scintillations  would  allow  one  to  design 
methods  to  compensate  for  them.  A  possible  long-range  goal 
would  be  the  ability  to  predict  when  and  under  what  physical 
conditions  that  scintillations  occur.  Evidently  this  is  being 
done  now,  but  only  very  roughly  (Fremouw  and  Wittwer,  1984:100) 
The  study  of  ionospheric  scintillations  is  of  interest  to  the 
civilian  community  as  well.  Applications  include  communica¬ 
tions  with  remote  sensing,  scientific,  and  commercial  satel¬ 
lites  in  high  inclination,  i.e.,  polar,  orbits;  this  includes 
the  recently  activated  search-and-rescue  satellites.  For  this 
and  other  reasons,  the  National  Research  Council  of  Canada  has 
agreed  to  the  collection  and  sharing  of  HILAT  data  with  the 
United  States  (Fremouw,  1983:6). 

As  was  earlier  mentioned,  it  has  been  proposed  that  low- 
energy  electron  precipitation  in  the  auroral  oval  is  the  pri¬ 
mary  source  of  plasma  structure  in  the  ionosphere  (the  F  region 
to  be  more  specific).  Basu  et_  al_.  performed  a  case  study  of 
the  formation  of  one  of  these  so  called  "blobs"  by  low-energy 
(less  than  500  ev)  electron  precipitation  near  the  equatorward 
edge  of  the  auroral  oval  during  an  auroral  substorm.  The  pre¬ 
cipitation  event,  monitored  by  a  DMSP  satellite,  caused  an 


increase  in  the  background  ion  density  by  a  factor  of  2-3  at 
an  altitude  of  820  km,  bounded  by  regions  of  sharp  density 
gradients.  Coordinated  magnetometer  and  total  electron  count 
measurements  showed  that  the  event  lasted  approximately  ten 
minutes.  It  was  also  noted  that  the  number  of  medium  scale 
(1  km  or  larger)  irregularities  responsible  for  phase  and 
amplitude  scintillations  increased  as  well.  In  this  study, 
scintillation  data  was  obtained  from  244  MHz  transmissions 


from  a  Fleetsatcom  satellite  in  geostationary  orbit  (Basu 
et  al . ,  1983:1151-1152)  . 

Since  no  electric  field  (or  field-aligned  current)  measure 
ments  were  available,  it  was  difficult  to  ascertain  the  genera 
tion  mechanisms  of  the  scintillation-producing  irregularities. 
However,  since  the  density  enhancement  was  embedded  in  the 
diffuse  aurora  region  with  known  field-aligned  Birkeland  cur¬ 
rents,  there  remains  the  possibility  that  these  same  Birkeland 
currents  are  destabilizing  elements  leading  to  irregularity 
formation  (Basu  et^  al_.  ,  1983:1163). 

Also,  the  data  indicated  the  possibility  of  a  fluctuating 
zonal  electric  field  with  alternating  eastward  and  westward 
components.  Since  the  medium  scale  irregularities  have  suf¬ 
ficiently  long  lifetimes,  a  fluctuating  electric  field  could 
introduce  instability  as  well.  There  is  also  the  possibility 
that  irregularities  may  grow  out  of  structured  low-energy  par¬ 
ticle  precipitation  and  irregular  field-aligned  currents 
(Basu  et_  al .  ,  1983:1163). 


With  the  advent  of  HILAT,  with  its  driftmeter,  magnetom¬ 
eter,  and  higher  spatial  resolution  of  particle  measurements, 
hopefully  more  can  be  learned  about  scintillation-producing 
irregularities  associated  with  F  region  density  enhancements 
(Basu  et_  al_.  ,  1983:1163). 

Problem  Statement 

The  main  thrust  of  this  thesis  effort  is  an  extension  of 
the  study  done  at  AFGL.  An  important  question  to  be  answered 
is,  "How  much  of  a  contribution  do  low-energy  electron  precip¬ 
itation  events  make  to  the  observed  scintillation  levels  occur¬ 
ring  in  the  auroral  zone?"  By  performing  a  survey  of  HILAT 
passes  at  different  magnetic  local  times,  one  could  determine 
how  often  these  enhancements  occur,  whether  or  not  there  is 
any  local  time  dependence,  and  how  often  their  observation  is 
accompanied  by  scintillations.  In  addition,  magnetic  field 
characteristics  associated  with  these  enhancements  could  be 
examined . 

Ob  i  ect ives  of  the  Research 

The  first  objective  is  to  investigate  the  above  mentioned 
possibility  of  magnetic  local  time  dependence  of  enhancement 
occurrence.  Second,  any  magnetic  field  characteristics  asso¬ 
ciated  with  these  enhancements  will  be  studied. 

Third,  the  electron  peak  number  densities  associated  with 
the  observed  enhancements  will  be  established.  Once  this  is 
done,  the  distribution  of  these  number  densities  could  be 
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plotted,  and  an  average  peak  number  density  value  might  then 
be  available. 


As  the  fourth  objective,  an  attempt  to  correlate  amplitude 
and  phase  scintillations  described  by  measured  values  of  the 
S4  and  PRMS  indices  with  the  position  and  magnitude  of  the 
particle  enhancement  will  be  made. 

Scope 

This  analysis  of  the  relationship  between  electron  enhance 
ments  and  scintillations  is  limited  to  electron  energies  and 
number  densities  measured  by  the  HILAT  electron  spectrometer 
(20ev-20kev) .  Radio  scintillation  S4  and  PRMS  data  is  from 
the  HILAT  multifrequency  radio  beacon;  the  frequencies  to  be 
used  are  138  MHz  and  413  MHz.  Magnetometer  data  will  also  be 
from  the  HILAT  vehicle  as  well.  The  time  interval  of  interest 
in  this  effort  covers  approximately  the  period  beginning  in 
early  January  1984  up  to  and  including  mid-March  1984. 


II.  Methodology  and  HILAT  Data 


Methodology 

By  using  expanded  electron  spectrometer  (herein  known  as 
the  J-sensor)  plots,  the  number  and  location  in  terms  of  mag¬ 
netic  local  time  (MLT)  of  the  enhancements  can  be  determined. 
Using  24  one-hour  "bins",  the  distribution  of  the  number  of 
enhancements  versus  MLT  can  be  obtained. 

The  remaining  analysis  will  use  only  a  selected  number  of 
representative  enhancements.  For  each  of  these  enhancements, 
any  corresponding  changes  in  the  magnetic  field  (derived  from 
expanded  magnetometer  data)  will  be  noted.  Using  printouts 
from  an  interactive  computer  program  from  AFGL,  the  peak  inten¬ 
sity  of  the  enhancements  should  be  easily  available.  With  this 
information,  a  distribution  of  peak  number  densities  can  be 
plotted.  Measuring  the  width  of  the  enhancements  in  terms  of 
time  (from  J-sensor  plots)  will  give  the  physical  size  of 
these  "blobs",  since  the  satellite's  orbital  velocity  is  known. 

The  AFGL  interactive  program  can  render  individual  number 
density  distribution  spectra  continuously  over  the  lifetime  of 
the  enhancement;  from  this  some  possible  insight  on  enhance¬ 
ment  formation  and  termination  might  be  gained. 

Lastly,  computer  printouts  from  AFGL  containing  S4  and 
PRMS  data  will  be  used  to  attempt  correlation  of  these  two 
parameters  with  enhancement  locations  in  MLT  using  available 
statistical  software  packages. 


The  data  preparation,  analysis,  and  results  for  each  sub¬ 
objective  can  be  found  in  Chapter  IV  (Graphical  Distribu¬ 
tion  of  Enhancements),  Chapter  V  (Magnetometer  Plots),  Chap¬ 
ter  VI  (Graphical  Distribution  of  Peak  Intensities  and 
Descriptive  Statistics),  and  Chapter  VII  (Discriminant  Analy¬ 
sis  of  Enhancements  and  Their  Relationship  to  Scintillations). 

Data 

Spacecraft  Descript  ion  and  Orbit .  The  HILAT  satellite, 
also  known  as  P83-1,  is  a  modified  Navy  Transit  satellite  modi 
fied  by  the  John  Hopkins  Applied  Physics  Laboratory  in  order 
to  perform  scientific  work.  The  spacecraft  was  launched  on 
27  June  1983  from  Vandenberg  Air  Force  Base  by  a  Scout  rocket; 
the  mission  orbit  has  a  perigee  of  800  km  and  an  apogee  of 
807  km  (Potocki,  1984:107). 

It  was  necessary  to  have  an  inclination  that  would  result 
in  a  high  precession  rate  (in  order  to  sample  all  possible 
local  times)  and  good  coverage  of  the  magnetic  polar  region. 
Consequently  the  satellite  was  placed  in  orbit  with  an  82 
degree  inclination.  This  inclination  also  allows  for  direct 
comparison  of  data  obtained  from  the  DMSP  spacecraft  as  well, 
since  the  two  vehicles  would  travel  in  coplanar  orbits  during 
certain  parts  of  the  year  (Fremouw,  1983:19).  HILAT's  orbit 
has  a  period  of  101.95  minutes  with  the  orbital  plane  pro¬ 
cessing  at  a  rate  of  -7.5  minutes  per  day,  or  24  hours  in  a 
period  of  six  months.  The  orbital  velocity  is  7.436  km/sec. 
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with  an  overhead  pass  requiring  IS. 7  minutes  to  travel  from 
horizon  to  horizon  (Fremouw,  1983:19). 

HILAT  has  no  on-board  data  storage  capability;  all  science 
and  telemetry  data  is  transmitted  to  receiving  stations  in 
real  time.  Three  of  the  stations  (Tromso,  Norway,  and 
Ft.  Churchill,  Canada,  both  of  which  are  in  the  auroral  zone, 
and  Sondre  Stromfjord,  Greenland,  which  is  near  the  polar  cusp) 
are  fixed,  and  the  third,  located  near  Seattle,  Washington, 
is  transportable  (Fremouw,  1983:20). 

Multifrequency  Radio  Beacon .  This  HILAT  experiment  is  in 
actuality  the  modified  spare  payload  from  the  Wideband  mission 
mentioned  in  the  "Background"  section.  Its  transmission  fre¬ 
quencies  are  137.676  MHz,  378.609  MHz,  413.028  MHz,  447.447 
MHz,  and  1239.084  MHz  (L-band)  (Fremouw,  1983:12). 

The  first  frequency  lies  in  the  very  high  frequency  range 
(VHF)  and  the  next  three  lie  in  the  ultrahigh  frequency  (UHF) 
range.  All  four  are  used  for  complex-signal  scintillation 
measurements,  whereas  the  L-band  frequency  is  used  only  for 
amplitude  scintillation  measurements.  The  L-band  frequency 
also  serves  as  a  phase  reference  for  the  VHF  and  UHF  signals, 
as  well  as  the  telemetering  signal  to  the  ground  stations 
(Fremouw,  1983:12). 

Although  not  explicitly  used  in  this  analysis,  the  three 
UHF  signals  are  also  used  to  obtain  total  electron  count  (TEC). 
The  change  in  TEC  is  proportional  to  the  measured  difference 
in  phase;  the  phase  difference  is  in  turn  proportional  to 
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frequency  and  changes  in  the  phase  path.  The  changes  in  phase 
path  can  be  due  to  either  a  change  in  the  real  path  (where 
the  change  in  phase  would  be  proportional  to  frequency)  or 
small  changes  in  the  refractive  index  (where  the  change  in 
phase  is  inversely  proportional  to  frequency).  By  making 
measurements  on  adjacent  UHF  frequencies,  it  will  be  possible 
to  tell  if  the  phenomenon  is  due  to  changes  in  height  or  a 
change  in  TEC  (Davies,  1969:229-230). 

The  L-band  and  the  UHF  transmissions  have  left  circular 
polarization.  This  is  done  to  minimize  interaction  between 
the  UHF  and  L-band  antennas,  which  are  nested  together  on 
the  earthward  surface  of  the  vehicle.  The  VHF  antenna  on 
the  other  hand  is  mounted  on  the  end  of  one  of  the  solar  pan¬ 
els.  The  effective  radiated  power  levels  near  the  center  of 
the  antenna  beams  are  given  as  26  dBm  at  VHF,  21  dBm  for  each 
of  the  UHF  signals,  and  30  dBm  at  L-band  (referenced  to 
1  milliwatt)  (Fremouw,  1983:14). 

J- Sensor .  The  J-sensor  represents  a  means  to  analyze  a 
plasma  irregularity  in  terms  of  measured  electron  flux.  This 
instrument  employs  six  cylindrical  curved  plate  electrostatic 
analyzers  arranged  in  three  pairs.  Each  analyzer  consists  of 
three  components:  an  aperture  entrance,  a  set  of  two  con¬ 
centric  cyl  indr ical ly  curved  plates,  and  a  set  of  channeltron 
detectors.  A  channeltron  is  an  electron  multiplier  consist¬ 
ing  of  a  hollow  lead  glass  tube  with  a  large  potential  dif¬ 
ference  between  the  two  ends;  electrons  impacting  on  the  front 


end  eject  secondary  electrons  which,  after  acceleration  by 
the  applied  electric  field,  strike  the  walls  of  the  tube, 
thereby  ejecting  even  more  secondary  electrons.  This  cascade 
of  electrons  renders  a  measurable  pulse  at  the  channeltron 
output  which  is  then  counted  by  the  detector's  electronics 
(Hardy  e_t  al  .  ,  1984:125,130). 

The  cyl indr ical ly  curved  plates  of  the  analyzer 
selects  only  those  electrons  with  a  given  energy  corresponding 
to  an  applied  voltage  across  them.  If  the  velocity  of  the 
incoming  electron  and  the  applied  electric  field  force  are 
matched,  the  electron  will  follow  a  trajectory  that  causes  it 
to  impact  the  front  end  of  the  channeltron  (Hardy  trt  al . , 
1984:125) . 

For  each  analyzer,  one  set  of  cylindrical  plates  has  a 
60  degree  radius  of  curvature,  and  the  other  set  has  a  127 
degree  radius  of  curvature.  The  latter  detector  measures 
electrons  in  eight  semi logr ithmical ly  spaced  channels  between 
20  and  632  ev  with  an  acceptance  angle  of  4  by  6  degrees. 

The  former  measures  electrons  in  a  likewise  manner,  but  with 
an  energy  range  of  632  to  20,000  ev  and  an  acceptance  angle 
of  2  by  9  degrees  (Hardy  et_  al .  ,  1984:125). 

The  high-energy  detector  has  a  geometric  factor  at  the 
peak  of  its  response  curve  of  .001  square  cm-steradian  with  a 
fractional  energy-spectral  channel  width  of  9  percent.  For 
the  low-energy  detector,  the  geometric  factor  is  .0004  square 
cm-steradian  with  a  corresponding  channel  width  of  13  percent 


(Hardy  et_  aK  ,  1984:126). 

The  three  analyzer  pairs  are  oriented  so  that  one  pair 
measures  electrons  incident  from  the  satellite's  local  zenith 
(parallel  with  the  magnetic  field  lines  at  high  latitudes). 

The  second  pair  is  angled  40  degrees  from  local  zenith,  and 
the  third  pair  measures  electrons  from  the  local  nadir  (up- 
welling  energetic  electrons)  (Fremouw,  1983:16). 

The  J-sensor  was  also  designed  to  work  in  three  different 
modes.  In  mode  1,  only  the  low-energy  zenith  analyzers  are 
used.  Twenty-four  8-point  spectra  over  the  low  energy  range 
are  gathered  in  one  second;  this  results  in  a  spatial  resolu¬ 
tion  of  300  meters  (that  is,  the  satellite  travels  300  m  in 
the  time  necessary  for  one  spectrum).  In  mode  2,  all  channels 
of  both  zenith  sensors  are  sampled,  giving  12  16-point  spectra 
per  second,  and  a  spatial  resolution  of  600  meters.  Mode  3 
operation  uses  all  channels  in  all  three  pairs;  this  gives 
four  16-point  spectra  per  second  in  all  three  directions, 
and  a  spatial  resolution  of  1.8  km  (Hardy  et  arl.,  1984:126). 

Magnetometer .  The  three-axis  fluxgate  magnetometer  was 
originally  part  of  the  attitude  determination  system  of  the 
Transit  spacecraft.  By  adding  a  microprocessor,  improved 
resolution  and  sampling  rates  became  possible  (Potemra  et_  al .  , 
1984  : 120) . 

The  three  magnetometer  sensors  are  arranged  so  that  two 
are  on  two  different  solar  panels,  with  the  remaining  sensor 
affixed  to  the  spacecraft  body.  If  we  define  a  right-handed 


satellite  centered  coordinate  system  with  positive  x  in  the 
direction  of  the  velocity  vector  and  positive  z  pointing  down¬ 
ward,  the  solar  panel  sensors  are  in  the  positive  x  and  y 
directions.  At  the  polar  regions,  the  resulting  x-y  plane 
is  almost  perpendicular  to  the  magnetic  field  lines;  conse¬ 
quently,  these  sensors  will  detect  magnetic  disturbances 
caused  by  field-aligned  currents.  The  x  and  y  sensors,  when 
aligned  in  this  fashion,  also  minimize  the  effect  of  magnetic 
fields  associated  with  currents  flowing  in  the  solar  panels. 

The  z  sensor  will  essentially  measure  the  strength  of  the  total 
magnetic  field  when  over  the  pole,  since  the  satellite's  z 
direction  will  be  almost  parallel  to  the  field  lines  (Potemra 
et  al. ,  1984:120) . 

The  HILAT  magnetometer  provides  a  magnetic  field  resolu¬ 
tion  of  13.4  nanoteslas.  The  magnetometer  processor  can  pro¬ 
vide  20  magnetic  field  samples  per  second;  as  a  result,  field- 
aligned  currents  can  be  detected  with  a  spatial  resolution  of 
400  meters  (Potemra  et  al . ,  1984:121). 
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III.  Literature  Review 

Radio  Frequency  Scintillation 

One  publication  discovered  early  in  the  literature  review 
was  a  compilation  of  papers  edited  by  Aarons;  although  many 
of  these  papers  discussed  scintillations  as  they  relate  to 
distant  radio  stars,  many  of  the  concepts  and  techniques  were 
easily  applicable  to  radio  transmitting  satellites.  Except 
for  the  paper  by  Briggs  and  Parkin,  all  works  referenced  in 
this  section  are  from  Aarons'  book. 

A  paper  by  Harrower  was  concerned  with  the  theoretical 
interpretation  of  radio  star  scintillations.  The  general 
characteristics  of  scintillations  and  the  associated  iono¬ 
spheric  irregularities  were  discussed,  along  with  methods  in 
which  the  size,  shape,  orientation,  and  drift  velocity  of  the 
irregularities  causing  them  could  be  found  (Aarons,  1963:38-64). 

Aarons  himself  presented  results  done  with  low  elevation 
angle  scintillations  obtained  through  observations  made  at  the 
Sagamore  Hill  Radio  Observatory.  One  of  the  discoveries  indi¬ 
cated  that  scintillations  are  produced  both  in  the  troposphere 
and  the  ionosphere  (Aarons,  1963:68-75). 

Scintillations  of  tropospheric  origin  were  observed  to 
have  long  periods  (20  sec-2  min)  with  good  correlation  of 
individual  fluctuations  measured  at  different  frequencies. 

Since  the  radio  source  (Cygnus),  observed  in  the  50  to  400  MHz 
range,  rose  in  the  sky  with  its  line  of  sight  in  the  F  region 
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passing  through  auroral  irregularities,  scintillations  were 
always  evident  (Aarons,  1963:75-79). 

The  correlation  of  scintillations  with  geomagnetic  activ¬ 
ity  was  investigated  as  well.  On  quiet  days,  scintillations 
associated  with  the  rising  of  Cygnus  had  the  same  power  den¬ 
sity  spectra  (a  semilogrithmic  plot  of  received  power  versus 
time)  among  the  observing  frequencies,  but  the  recorded  peaks 
and  valleys  were  not  always  correlated  in  time.  Differential 
refraction  in  the  50-200  MHz  range  was  believed  responsible 
for  this  effect.  A  single  scattering  process  (thin  layer  of 
irregularities)  was  thought  to  be  a  possibility  as  well 
(Aarons ,  1963:79). 

On  geomagnetically  active  days,  multiple  scattering  (dif¬ 
fraction  process)  caused  the  power  spectra  for  the  different 
observing  frequencies  to  differ.  Here  the  scale  of  the  ground 
projection  of  the  observed  irregularities  was  such  that  a  low 
scintillation  index  at  54  M'lz  and  a  high  index  at  244  MHz 
resulted,  an  example  of  the  inversion  effect  (Aarons,  1963: 
79-81) . 

This  effect,  which  was  discussed  by  Aarons  et_  aT.  in 
another  publication,  manifests  itself  when  the  correlation 
distance  on  the  ground  under  multiple  scattering  conditions  is 
smaller  than  the  antenna  diameter.  Several  peaks  and  nulls 
exist  across  the  antenna  dish,  and  due  to  a  lack  of  coherence, 
sc int  i  1  lat  io*.  peaks  and  nulls  combine  to  produce  a  composite 
low  index.  As  the  index  of  refraction  is  inversely  proportional 


• '  j>-^jr»>TygT.’%-g^A  '.v'.v'.v  v,«r,TO  ^  ».’  H"H5V.\g  WWVT.I  tv  »r  *v  jrxru  r. 


rsprrawr*j»^*^  .• 


to  frequency  squared,  higher  frequencies  experience  less  scat¬ 
tering;  therefore,  the  correlation  distances  on  the  ground  are 
almost  equal  to  the  actual  size  of  the  irregularities.  As  a 
result,  the  scintillation  index  appears  to  increase  with  fre¬ 
quency,  which,  as  will  soon  be  shown,  is  usually  not  the  case 
(Aarons,  1963:81-83). 

In  another  paper  by  Mass  where  various  satellite  techniques 
for  studying  propagation  in  the  ionosphere  were  surveyed, 
scintillations  were  considered  as  one  of  the  major  disruptive 
effects.  According  to  Mass,  amplitude  scintillations  are  the 
most  detrimental  with  respect  to  communications.  It  was  also 
pointed  out  that  good  correlation  exists  between  the  occurrence 
of  satellite  amplitude  scintillations  and  amplitude  scintil¬ 
lations  associated  with  radio  stars  (Aarons,  1963:283). 

Frihagen  presented  findings  on  the  altitude  of  the  scin¬ 
tillation  producing  irregularities;  his  results  indicated  that 
this  region  was  between  300  and  550  km.  The  size  and  the 
orientation  were  also  determined,  with  the  irregularities  in 
this  case  being  at  least  1.3  km  long  and  no  more  than  3.6  km 
wide.  They  were  also  anisotropic,  i.e.,  they  were  aligned 
with  the  magnetic  field  lines  (Aarons,  1963:342). 

The  final  paper  reviewed  from  Aarons'  book  was  by  Lizka. 

His  results  taken  from  various  satellite  studies  of  the  auroral 
zone  indicated  that  auroral  sporadic-E  regions  produced  scin¬ 
tillations.  Also,  observations  from  Explorer  7  showed  that 
the  frequency  of  scintillation  occurrence  decreased  as  distance 
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from  the  auroral  zone  increased  (Aarons,  1963:345-346). 

After  reviewing  the  previously  cited  works,  it  became  ap¬ 
parent  that  no  common  scintillation  index  was  in  use  (at  least 
in  the  early  1960s,  when  most  of  the  papers  were  written). 

Since  the  S4  amplitude  index  and  the  PRMS  are  the  scintilla¬ 
tion  parameters  of  interest  in  this  thesis  effort,  it  was 
necessary  to  find  some  source  in  the  literature  that  would 
explain  just  what  exactly  is  being  measured  by  these  parameters. 

A  study  on  the  variation  of  radio  star  and  satellite 
scintillations  with  zenith  angle  by  Briggs  and  Parkin  provided 
actual  mathematical  expressions  for  both  PRMS  and  the  S4 
index.  They  start  by  considering  fluctuations  on  a  wave  as  it 
passes  through  an  irregular  medium;  in  this  instance,  the 
medium  is  considered  as  a  thin  diffracting  screen.  Since  the 
frequencies  most  likely  to  be  used  will  suffer  little  iono¬ 
spheric  absorption,  the  screen  will  produce  across  the  emerg¬ 
ing  wavefront  phase  variations  only  and  no  amplitude  varia¬ 
tions.  Only  when  the  wave  propagates  beyond  the  screen  do 
amplitude  fluctuations  develop  (Briggs  and  Parkin,  1963: 

339-340) . 

After  discussing  the  geometry  of  the  problem  (Fig.  1)  and 
the  propagation  of  a  wave  through  the  medium  (see  Fig.  2  for 
a  schematic  representation  of  an  irregularity),  the  general 
equation  for  the  root  mean  square  fluctuation  of  phase  pro¬ 
duced  in  a  wave  traveling  along  the  z-direction  is  (Briggs 
and  Parkin,  1963:345) 


FIG  2. 


Diagram  of  an  Irregularity 
(Briggs  and  Parkin,  1963:  343) 


.  PRMS^TT'^X  UN:^*r#secxH*xsi^4'  +  cos1  S')  -{1) 

where 

=  classical  electron  radius 
X  =  wavelength  of  signal 
*N  =  mean  excess  electron  density 
r0  =  traverse  radius  of  irregularity 
e*  =  ratio  of  irregularity  length  to  width 
The  variable  *  h  is  the  thickness  of  the  irregularity,  and  the 
angle  i  is  the  angle  of  incidence  that  the  radio  beam  makes 
with  respect  to  the  irregularity's  surface  (Fig.  1).  The 
angle  T  appears  in  Fig.  2;  it  is  the  angle  between  the  mag¬ 
netic  field  vector  H0and  the  direction  of  wave  propagation. 

The  only  factors  that  vary  with  zenith  angle  ©  are  i  and 
'P  .  The  expression  is  simplified  by  assuming  that  fr ,  h, 
and  r0  are  independent  of  O  : 

PRMS  =  K  X  (sec  iJVsin^t  Cos'*  S')  *  (2) 

where  K  is  a  constant  of  proportionality.  For  isotropic  (non- 
aligned  with  the  magnetic  field)  irregularities,  further  sim¬ 
plification  shows  PRMS  varying  asX(sec  i)^  (Briggs  and  Parkin, 
1963:345)  . 

Briggs  and  Parkin  then  reference  a  paper  by  Mercier  to 

explain  how  the  amplitude  scintillation  index  S  behaves. 

2 

S  is  used  to  measure  scintillation  depth,  and  it  is  expressed 


s2  =  <R*>a(<R*>  -  <«.*>*) 


‘  *  *  *  ’  >  *  »*•  -*•  k‘»  A  k'i  A  A  _ V  -‘.V*  ■  *  s."  \ 


(3) 


where  R  is  wave  amplitude.  The  quantity  S  also  varies  as  a 
function  of  X  z/r^  (Fig.  3);  note  that  near  the  origin  S  is 
proportional  to  z.  Briggs  and  Parkin  refer  to  this  region 
as  the  "near  zone".  Further  away  from  the  screen  S  approaches 
a  limiting  value,  and  S  becomes  independent  of  z;  this  is  the 
"far  zone".  The  limiting  value  of  S  far  from  the  screen  is 

S  =  {l-exp(-2 (PRMS) 2 } ^  (4) 

as  determined  by  Mercier  (Briggs  and  Parkin,  1963:346). 

For  small  values  of  PRMS,  S  can  be  expressed  as 

s  =  rr  (prms)  (i  +  ttV/h-xV)’7*  (5) 

where  the  limiting  value  far  from  the  screen  is  (PRMS) 
(Briggs  and  Parkin,  1963:346-347). 

Mercier's  paper  applied  to  isotropic  irregularities.  For 
the  more  general  anisotropic  case,  James  put  forth  the  follow¬ 
ing  relationship 

S  =  iT(PRMS)  { 1- (cosu1cosu2)  ^cos(l/2)(u^+U2) }  ^  (6) 

where  tan  Uj  =  2Xz/TT  r^  and  tan  U2  =  2>z/TY  P*  r*  .  The  param¬ 
eter  &  is  an  axial  ratio  of  the  phase  pattern  expressed  as 
(Briggs  and  Parkin,  1963:344,347) 

Y  f  Cosa  (7) 

For  the  case  of  scintillations  of  satellite  transmissions 
caused  by  anisotropic  irregularities,  the  PRMS  value  can  be 


found  with  Eq  (2),  and  B  can  be  determined  from  Eq  (7).  Th< 
scintillation  depth  can  be  determined  from  either  Eq  (6)  or 
Fig.  4;  however,  z  in  Eq  (6)  must  be  replaced  by 


{z^  +  l2)  *  The  values  of  z^  and  z^  (Fig.  1)  are  deter¬ 
mined  by  (Briggs  and  Parkin,  1963:358) 

2,=  (itfcos*©  +  U  .(8) 

R*  COSa©  +  -  -Z,  .  (9) 

In  the  case  where  two  frequencies  are  being  used  by  the 
satellite,  the  ratio  of  the  scintillation  depth  is  indepen¬ 
dent  of  PRMS.  This  ratio  is  given  as 

-V±  >/a 

SU.J/SU^’  X./X^C  1  +  tfr*/4X^£).  CIO) 

where  Z  =  z1z2/ (Zj  +  z^)  (Briggs  and  Parkin,  1963:360). 

The  authors  go  on  to  say  that  for  small  values  of  r  in 

o 

the  far  zone,  the  limiting  case  of  Eq  (10)  is  (Briggs  and 
Parkin,  1963:355) 

S(X,)/S(X>)  =  X(/X^  (11) 

If  rQ  is  large  (observations  made  in  the  near  zone),  we  have 

S(Xl)/S(Xa)  =  (X./xj2  (12) 

As  was  mentioned  previously,  when  this  paper  was  written 
different  measures  of  scintillation  depth  were  in  use. 

Briggs  and  Parkin  felt  that  their  theoretical  calculations 


JV  \  'w  -  '  "V  --1^. 


would  not  be  comparable  to  experimental  data  unless  some  rela¬ 
tionship  between  the  two  were  specified  (Briggs  and  Parkin, 
1963:362) . 

When  radio  waves  are  received  from  a  satellite,  recorder 
deflection  is  proportional  to  the  signal  amplitude.  Two 
measurements,  the  mean  deviation  (SI)  and  the  root -mean-square 
deviation  (S 2 )  are  possible.  SI  and  S2  can  be  expressed  as 
(Briggs  and  Parkin,  1963:362) 


SI  =  <R») 

(13) 

S2  =  -  <R»;1},/3' 

(14) 

When  observing  radio  star  scintillations,  recorder  deflec 

2 

tion  is  proportional  to  R  ,  or  received  power.  The  two 
measurements,  S3  and  S4,  are  expressed  as  (Briggs  and  Parkin, 
1963:363) 


S3  = 

<**>*(<  Ra- 

<fO*» 

(15) 

_  •*  i  , 

S4  = 

<  RJ>  {<* - 

<  p?» } 

(16) 

Although  not 

readily  apparent,  Eq  (16) 

and 

Eq  (3)  give 

similar 

results , 

i .e . ,  S=S4 . 

The  authors  were 

then 

able  to 

develop 

empirical 

re lat ionship 

is  that  relate 

SI, 

S2  , 

and  S3 

to  S4  (Briggs  and 

Parkin,  1962 

1:364-365) 

SI  = 

0.42  S4 

(17) 

S2  = 

0.52  S4 

(18) 

S3  =  0.73  S4 


(19) 


■*.  •  X. 


Based  on  private  communications  with  Dr.  Sunanda  Basu, 

PRMS  is  inversely  proportional  to  frequency  for  all  frequen¬ 
cies,  in  agreement  with  Eq  (1).  This  is  not  at  all  the  case 
with  S4,  as  can  be  seen  in  either  Eq  (5)  or  (6).  A  good  rule 
of  thumb  according  to  Dr.  Basu  is  that  S4  is  inversely  pro¬ 
portional  to  frequency  to  the  three-halves  power,  and  only 
for  S4  values  less  than  approximately  .5  or  .6  (Basu,  1984: 
personal  interview). 

Sources  of  Plasma  Structure 

There  are  many  mechanisms  believed  responsible  for  the 
creation  of  plasma  irregularities;  in  addition  to  precipitat¬ 
ing  electrons  what  follows  is  a  sample  of  some  of  the  studies 
done  and  some  of  the  resulting  theories  concerning  irregular¬ 
ity  formation. 

One  of  the  earlier  explanations  was  offered  by  Dagg  con¬ 
cerning  current  flow  from  an  induced  electromagnetic  force. 

At  the  height  of  the  D  region,  ions  and  electrons  move  with 
neutral  particles  transversely  to  the  goemagnetic  field  lines. 
This  motion  is  transferred  along  the  field  lines  to  the  higher 
ionospheric  layers,  which  is  between  100  to  125  km  high. 
Pressure  variations  in  the  lower  atmosphere  effect  electrons 
and  ions  in  the  E  layer.  The  conductivity  of  the  ionosphere 
at  this  altitude  is  the  same  both  parallel  and  transverse  to 
the  field  lines;  consequently,  the  resulting  charged  particle 
motion  mimics  a  moving  conductor  in  a  reciprocating  dynamo. 
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The  resulting  current  flow  manifests  itself  as  a  change  in 
the  magnetic  elements  (Aarons,  1963:40). 


TIP .  Many  researchers  have  investigated  traveling  iono¬ 
spheric  disturbances  (TID)  ,  which  are  essentially  large  elec¬ 
tron  density  waves  (more  than  1500  km  wide)  propagating  in 
the  F  region.  The  origin  of  these  waves  may  be  linked  to 
acoustic  gravity  waves  propagating  in  the  thermosphere  (200  to 
500  km) .  Internal  gravity  waves  are  atmospheric  disturbances 
with  periods  ranging  from  10  minutes  to  several  hours.  A 
wave  duct  may  result  between  a  reflecting  region  and  the  sur¬ 
face  where  these  gravity  waves  are  reflected  in  a  thermally 
structured  atmosphere.  In  cases  where  reflection  is  not  quite 
total  but  still  strong,  some  wave  energy  can  "leak"  into  the 
higher  ionosphere  (Stevens,  1973: (6-1) - (6-3) )  .  Friedman, 
Elkins  and  Slack,  and  Yeh  have  contributed  papers  concerning 
the  TID  phenomenon  and  its  relationship  to  plasma  structures 
(Friedman,  1966:1033-1054;  Elkins  and  Slack,  1969:421-439; 


Yeh,  1972:709-719). 

Sporadic  jj.  Under  certain  conditions,  irregularities  are 
aligned  with  the  geomagnetic  field  which  occur  in  the  E  region 
near  the  magnetic  equator.  A  mechanism  put  forth  by  Farley 
responsible  for  these  irregularities  is  the  two-stream,  or 
plasma  ion  wave  instability.  A  plasma  with  two  or  more  inter¬ 
connecting  streams  o ;  charged  particles  will  be  unstable  if 


the  mean  velocity  of  the  particles  of  one  stream  is  much 
greater  than  the  mean  velocity  of  the  particles  of  the  other 


(Stevens,  1973 : (6- 9) ) . 


Farley  put  forth  the  arguemnt  that  irregularities  should 
occur  in  the  equatorial  electrojet  because  of  two-stream  in¬ 
stabilities.  Farley  also  argued  that  similar  irregularities 
appearing  in  the  polar  regions  are  caused  by  auroral  electro¬ 
jets  (Farley,  1963:6083-6097).  Others  that  have  conducted 
sporadic  E  research  are  Ireland,  Preddey  (Ireland  and  Preddey, 
1967:137-148)  and  Goodman  (Goodman,  1967:607-612). 

Another  instability  in  the  equatorial  electrojet  is  the 
gradient  drift,  or  E  X  B  instability.  Here  E  is  the  ambient 
background  electric  field  and  B  is  the  geomagnetic  field. 
Farley  and  Balsey  studied  the  dependency  of  electrojet  irreg¬ 
ularities  on  magnetic  field  and  an  associated  critical  value 
(Farley  and  Balsey,  1973:2  2  7-239).  Sudan  et_  a_l .  have  studied 
this  instability's  effect  as  well  on  small  scale  irregulari¬ 
ties  as  observed  by  radar  (Sudan,  1973:240-248). 

Spread  £ .  This  phenomenon  of  the  F  region  is  character¬ 
ized  as  an  apparent  thickening  found  in  the  reflecting  re¬ 
gions  believed  to  be  caused  by  vertical  electron  density 
inhomogen  iet  ies .  Although  by  definition  spread  F  is  observed 
below  the  critical  frequency  of  the  F  region,  many  feel  these 
irregularities  are  responsible  for  scintillations  observed  at 
VHF  and  UHF  frequencies  (Stevens,  1973: (6-15)) . 

A  proposal  by  Liu  and  Yeh  employs  the  Rayleigh-Taylor 
instability  as  the  mechanism  which  causes  spread  F.  This 
instability  is  due  to  the  separation  between  a  heavy  upper 


fluid  and  a  lighter  lower  fluid,  which  is  unstable  in  a  gravi¬ 
tational  field.  An  analogous  situation  occurs  in  the  iono¬ 
sphere,  where  the  density  gradient  increases  with  altitude. 
After  the  application  of  a  perturbation,  the  system  will  be 
stable  at  the  lowest  possible  potential  energy,  with  extra 
energy  transformed  into  the  kinetic  energy  of  instability 
(Liu  and  Yeh,  1966:1407-1416).  Liu  and  Yeh  developed  a  dis¬ 
persion  relation  showing  that  when  the  height  gradient  of 
electron  density  is  positive  and  deep,  instabilities  may 
result.  Both  felt  experimental  data  supported  their  theory 
(Liu  and  Yeh,  1966:1283-1291).  Both  Farley  and  Reid  have 
independently  studied  the  spread  F  phenomenon  in  both  the 
polar  and  non-polar  regions  as  well  (Farley,  1960:869-877; 
Reid,  1968:1627-1640). 

Polar  Cap  VHF  Scintillation.  After  studying  40  MHz  scin¬ 
tillation  observations  taken  via  satellite  over  the  polar 
regions,  Frihagen  deduced  that  scintillation  depth  was  inde¬ 
pendent  of  geomagnetic  activity,  time  of  day,  or  satellite 
position.  He  went  on  to  develop  a  model  where  particle  precip 
itation  from  the  interplanetary  medium  resulted  in  F  region 
irregularities.  To  produce  the  observed  irregularities, 

11  _  O 

Frihagen  calculated  a  spatial  variation  value  of  10  m 
sec  ^  at  one  Kev  as  being  a  required  value  necessary  to  sus¬ 
tain  "blob"  production  (Frihagen,  1969:81-92). 


IV.  Electron  Data  Analysis 


Processing  Equipment  and  Software 

The  computer  used  to  perform  the  necessary  data  analysis 
was  the  Control  Data  Corporation  Cyber  6000,  which  is  pri¬ 
marily  used  by  the  Aeronautical  Systems  Division  but  is  also 
available  for  student  use  as  well.  The  primary  reason  why  the 
Cyber  was  selected  is  because  the  required  statistical  analy¬ 
sis  program  resides  in  the  Cyber  system. 

This  statistical  analysis  software  is  the  Statistical 
Package  for  the  Social  Sciences  (SPSS),  designed  by  Nie  et  al . 
originally  at  Standford  University  in  1965.  This  package 
was  selected  because  it  is  quite  user-friendly  and  it  is  well 
documented.  In  this  effort  the  two  subprograms  FREQUENCIES 
and  DISCRIMINANT  were  the  primary  SPSS  routines  used. 

Also,  enlarged  J-sensor  plots  and  distribution  spectra 
that  are  products  of  an  interactive  program  were  used  in  this 
research  effort.  This  program  was  designed  by  Dr.  Robert 
Burkhart  of  Emmanual  College,  which  is  closely  affiliated 
with  AFGL. 

Graphical  Distribut ion  of  Electron  Enhancements 

The  first  subobjective  of  the  research  is  to  see  how  the 
low-energy  electron  enhancements  are  distributed  in  magnetic 
local  time  (MLT) .  Since  this  and  the  remaining  subobjectives 
all  study  some  aspect  of  these  enhancements,  it  would  be 
prudent  at  this  point  to  define  just  what  constitutes  an 


enhancement  and  show  an  example. 

Before  doing  that,  we  need  to  define  the  three  measure¬ 
ments  EAVG,  JETOT,  and  JTOT  that  appear  on  a  J-sensor  plot. 

If  one  considers  the  differential  number  flux  J  (E . ) 

n  v  1 ' 

2 

(electrons/cm  -sec-ster-keV)  where  i  runs  from  1  to  16 

(16  channels)  and  E.^  is  the  central  energy  of  each  channel  in 

keV,  then  the  integral  number  flux  JTOT  is  expressed  as 

JTOT  =  Jn(E,HEi-El)tI.Jh(EiHCE;,1-E.1.1>/l)+J„(Ej(E;E;(20) 

'•‘X 

Similarly,  the  differential  energy  flux  is  JC(E.) 

E  1 

2 

(keV/cm  -ster-keV) ;  the  expression  for  the  integrated  energy 
flux  is 


Finally,  the  expression  for  the  average  energy  is  just  the 
integrated  energy  flux  divided  by  the  integrated  number  flux 
(JETOT/JTOT) .  On  the  J-sensor  plots  EAVG  is  expressed  in 
keV  (Hardy,  undated : 8- 9) . 

Fig.  5  is  an  example  of  the  J-sensor  plots  mentioned  in 
the  "Methodology”  section  in  Chapter  1.  This  data  is  from  the 
Tromso  site,  taken  from  1:43:47  to  1:53:45  U.T.  on  17  January 
1984.  Fig.  5  is  actually  an  expanded  portion  of  the  total 
pass;  the  time  of  interest  is  1:48:40  to  1:51:40.  Notice  the 
two  peaks  appearing  in  the  bottom  JTOT  (integrated  number  flux) 
panel  at  approximately  1:50:01  and  1:50:16.  Notice  also  that 
both  peaks  are  superimposed  on  what  is  known  as  the  diffuse 


aurora,  which  runs  from  approximately  1:49:50  to  about  1:50:55 
The  region  immediately  preceding  this  (1:49:10  to  1:49:50)  is 
the  so-called  discrete  aurora,  and  any  such  peaks  found  in 
this  region  are  not  considered  enhancements  as  defined  here. 

The  two  peaks  have  corresponding  troughs,  or  'cutouts'  in 
EAVG  (average  energy  in  Kev,  top  panel)  which  are  approxi¬ 
mately  an  order  of  magnitude  less  than  the  ambient  background 
energy.  In  summary,  a  low-energy  electron  enhancement  has  a 
JTOT  peak  in  the  diffuse  aurora  and  a  corresponding  drop  in 
EAVG. 

Distribution  spectra  (Fig.  7)  are  more  quantitative  indi¬ 
cators  of  what  constitutes  an  enhancement.  In  Fig.  6  we  have 
an  enhancement  recorded  at  Sondrestrom  on  21  January  1984  at 
6:52:20.  Fig.  7  shows  the  distribution  spectra  for  the 
entire  pass  (24731  sec  to  24803  sec  after  midnight  U.T.),  or 
converting  into  the  more  familiar  hour,  minutes,  seconds  for¬ 
mat,  6:52:11  to  6:53:23.  We  have  a  three-dimensional  plot  of 
the  logrithm  of  the  energy  in  ev  versus  the  logrithm  of  the 
differential  number  flux  Jn(E)  (electrons/square  cm/sec/ster/ 
ev)  versus  time.  Note  the  increase  in  the  flux  value  at  the 
lower  end  of  the  energy  axis  shortly  after  24738  seconds, 
corresponding  with  the  observed  enhancement.  This  is  far  more 
evident  if  we  examine  the  individual  distribution  spectra, 
which  appears  in  Fig.  8  and  Appendix  A. 

Since  we  are  working  with  Mode  3  data,  we  have  four  indi¬ 
vidual  distribution  spectra  per  second.  Looking  at  the  very 
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first  one  (Fig.  8)  starting  at  24739  seconds  we  have  a  bi- 
Maxwellian  distribution  curve  (actually  the  logrithm  of  this 
distribution  F(E)  versus  energy  in  ev) .  In  the  upper  left- 
hand  corner  we  have  the  parameters,  electron  density  n  and 
temperature  T,  that  characterize  both  the  low  energy  and  high 
energy  curve.  Immediately  below  that  we  have  for  each  J- 
sensor  channel  the  number  of  raw  counts  registered;  the  infor¬ 
mation  in  the  other  two  columns  were  not  used. 

The  quantity  AEGY  appearing  below  the  raw  counts  table  is 
the  average  energy  EAVG  in  Kev .  The  endpoints  for  the  two 
curves  are  specified  by  channel;  in  this  case  El  is  the  20  ev 
channel,  E2  (the  "breakpoint")  is  the  632  ev  channel,  and  E3 
is  the  20  Kev  channel.  By  inspecting  the  graph  we  see  that 
both  curves  are  fairly  Maxwellian  in  nature. 

Note  that  as  we  look  at  subsequent  individual  spectra, 
we  see  a  significant  increase  in  low  energy  counts  and  a  con¬ 
sequent  lowering  of  average  energy  EAVG.  At  the  enhancement's 
peak  we  have  the  maximum  count  in  the  low  energy  channels 
(20-385  ev) ,  with  a  minimum  value  of  EAVG  (1.04  Kev).  As  we 
approach  the  end  of  the  enhancement  we  get  a  substantial 
lowering  of  the  low  energy  counts  and  an  increasing  EAVG  value 
It  was  not  possible  to  get  distribution  spectra  for  every 
individual  enhancement,  so  only  a  few  selected  cases  appear  in 
Appendix  B. 

In  order  to  establish  the  MLT  distribution  of  the  enhance¬ 


ments,  passes  from  9  January  1984  to  13  March  1984  were 


examined  (not  on  consecutive  days;  see  Table  1).  A  total  of 
198  enhancements  were  obtained,  and  the  corresponding  hour  of 
corrected  magnetic  local  time  (also  obtained  from  the  J-sensor 
plots)  were  noted.  This  data  was  entered  into  the  Cyber  and 
an  SPSS  program  (FREQUENCIES)  was  executed. 

Basically,  FREQUENCIES  determines  frequency  distribution 
tables  and  descriptive  statistics.  In  this  first  subobjective 
we  are  interested  in  how  the  enhancements  are  distributed 
over  24  one-hour  "buckets”,  i.e.,  a  histogram.  (For  example, 
the  1100  MLT  "bucket"  means  the  time  interval  from  1100  MLT 
to  1159  MLT  inclusive).  This  is  also  easily  obtained  from 
FREQUENCIES  by  specifying  the  necessary  option. 

A  histogram  depicting  the  results  of  this  first  subobjec¬ 
tive  appear  in  Fig.  9.  Note  the  large  clustering  in  the 
0300  MLT  sector  (post -midnight ) ;  there  is  also  a  secondary 
peak  appearing  in  the  2300  sector  (pre-midnight). 

Notice  that  from  0700  to  1400  MLT  there  is  only  one 
enhancement  (0800)  .  The  results  indicate  that  at  least  for 
the  time  period  of  this  research,  where  all  MLT's  are  being 
uniformly  sampled,  enhancements  appear  in  the  nighttime  sec¬ 
tors,  with  a  sharp  dropoff  at  magnetic  dawn  and  a  sharp 
increase  shortly  before  magnetic  dusk. 


Table  1.  MLT  Distribution  Data 
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FIGURE  9.  Histogram  of  198  Electron  Enhancements 


V.  Magnetometer  Plots 


The  second  subobjective  of  this  research  effort  is  to  see 
if  the  appearance  of  an  electron  enhancement  causes  a  change 
in  the  geomagnetic  field  (or,  equivalently,  is  there  an  asso¬ 
ciated  Birkeland  current).  Time  did  not  permit  analysis  of 
all  198  enhancements,  so  seventeen  passes  consisting  of  73 
enhancements  were  selected  for  this  aspect  of  the  thesis. 

These  particular  passes  were  selected  in  such  a  way  to  get  as 
even  a  contribution  as  possible  from  each  month.  The  selected 

passes  also  had  two  or  more  enhancements;  due  to  time  limita¬ 

tions  it  was  important  to  use  as  many  as  possible  from  the 
few  passes  available  for  analysis. 

An  example  of  a  magnetometer  plot  appears  in  Fig.  10; 
what  one  would  expect  to  see  is  a  sharp  change  in  the  slope 
of  the  curve,  which  is  the  y-component  (east-west)  of  the  mag¬ 
netic  field  plotted  against  time,  that  is  coincident  with  an 

enhancement.  The  change  in  the  By  component  is  associated 
with  a  Birkeland  (magnetic  field-aligned)  current  by  the  fol¬ 
lowing  reasoning  borrowed  from  Potemra,  et_  al . 

The  Birkeland  current  J  produces  a  magnetic  field  in 
accordance  with  Ampere's  law, 

J=  (  M0)'  [curl  B]  (22) 

where  =  4  IT  x  10  henry/meter  and  B  is  the  geomagnetic 


field.  Since  magnetic  disturbances  associated  with  the 


MEASURED  ( MILL ICRUSS ) 


Birkeland  currents  are  in  the  B  direction  as  previously 
stated,  we  have 

J*  =  J»,  =  Kf'tdBy/dx]  (23) 

where  x  is  to  the  north,  y  to  the  east,  and  positive  z  is 
down  toward  the  earth  and  parallel  to  the  field  lines. 

Since  M.0 is  known,  and  since  the  HILAT  orbital  velocity 
dx/dt  =  8  km/sec,  J  can  be  written  as 

J„  =  (dx/dt)_1dBy/dt  (24) 

=  0.1  dBy/dt  pA/square  meter  (25) 

where  dBy/dt  is  in  nanoteslas  per  second.  The  magnetometer 
plot  shown  has  By  in  milligauss,  which  is  easily  convertible 
into  nanoteslas,  where  one  gauss  equals  10~4  teslas  (Potemra, 
1984:  121-122)  . 

There  is  one  problem  that  one  must  be  aware  of  in  working 
with  HILAT  magnetometer  data;  the  spacecraft  is  oscillating 
slowly  about  its  roll  axis  (less  than  a  degree  to  either  side) 
so  it  is  difficult  to  tell  if  the  slope  of  the  line  is  a 
spacecraft - induced  effect  or  if  it  is  a  genuine  magnetic  dis¬ 
turbance.  Since  we  are  looking  for  'cutouts'  in  the  curve 
that  correspond  to  an  enhancement  occurrence,  we  measure  the 
slope  of  the  cutout  with  respect  to  the  ambient  slope.  Fig.  10 
shows  an  example  of  what  we  seek  (Hardy,  1984 : personal  inter- 
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This  data  belongs  to  the  same  pass  measured  in  the  en¬ 
hancement  distribution  section  (17  January  1984,  1:43:47  to 
1:53:45  U.T.).  Notice  the  V-shaped  structure  (cutout)  appear¬ 
ing  between  6619  seconds  and  6622  seconds.  The  third  enhance¬ 
ment  of  this  pass  starts  at  6616  seconds  and  lasts  until  ap¬ 
proximately  6621  seconds.  In  this  case  the  ambient  slope  is 
essentially  zero  (from  6610  to  6625  seconds),  so  after  taking 
the  slope  of  the  negative,  or  left  half  of  the  V  (which  indi¬ 
cates  a  downward  current,  or  upward  flowing  electrons; 

Potemra,  1984:124),  we  get  a  current  density  of  13 . 3pA/square 
meter.  The  right  half  of  the  V  is  essentially  the  mirror  im¬ 
age  of  the  left  half;  the  current  density  has  the  same  magni¬ 
tude  but  opposite  direction  (upward  current,  downward  flowing 
electrons;  Potemra,  1984:124). 

In  the  example  pass  given,  the  first  two  enhancements  have 
no  corresponding  disturbance.  We  would  then  like  to  see  how 
many  enhancements  are  associated  with  magnetic  disturbances, 
and,  if  possible,  is  enhancement  size  proportional  to  the  size 
of  the  disturbance. 

Of  the  73  enhancements  studied  in  this  particular  subob¬ 
jective,  only  nine  had  corresponding  changes  in  By  with 
respect  to  time  (dB  /dt).  The  results  appear  in  Table  2;  here 
"Time"  is  the  beginning  time  of  the  enhancement,  and  "Time 
Duration"  is  the  measured  duration.  There  appears  to  be  some 
credence  to  a  possible  relationship  between  the  peak  number 
flux  of  the  low-energy  electrons  and  the  current  density 
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Table  2.  Magnetometer  Results 


Site 

(T/S) 

Date 

Time(U.T.) 

(HH: MM: SS) 

Time  Duration 

(sec) 

PEAK 

dBy/dt 

(nT/sec) 

J 

(  A/m2) 

g 

•i  -J.  '.J* 

T 

17-1 

01:50:16 

5 

58.4 

133 

13.3 

•  -  *  -  * 

T 

19-1 

17:38:24 

4 

1.97 

16.7 

1.67 

V  *-• 

T 

21-1 

01:46:12 

3 

1.75 

20.0 

2.00 

‘,V,4 

S 

23-1 

06:0026 

3 

14.8 

12.9 

1.29 

,v. 

T 

25-1 

20:04:34 

2 

38.8 

100 

10.0 

T 

9-3 

21:30:53 

5 

46.2 

68.2 

6.82 

T 

9-3 

21:31:47 

3 

3.6 

16.7 

1.67 

T 

9-3 

21:31:55 

2 

32.5 

3.25 

T 

9-3 

21:32:20 

4 

4.4 

5.3 

0.53 

.  v  J 

L"'-  **•  " 

(Date  is  in  "Day-Month"  format,  and  PEAK  is  in  units  of 
10®/cm2-sec-ster) 


T  =  Tromso 


S  =  Sondre 


However,  with  such  limited  statistics,  any  attempted  linear 
regression  would  yield  results  of  dubious  value.  The  nine 
points  appear  in  Fig.  11,  which  is  a  plot  of  PEAK  versus  J^. 
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Graphical  Distribution  of  Peak  Intensities 


and  Descriptive  Statistics 


In  this  third  subobjective,  we  are  interested  in  measuring 
the  peak  number  flux  (JTOT)  of  the  enhancements,  and  then 
studying  how  these  peak  values  are  distributed  along  with  some 
descriptive  statistics. 

For  this  subobjective,  the  same  seventeen  passes  with 
their  73  enhancements  were  used.  Enlarged  J-sensor  plots 
were  used  to  measure  the  peak  intensities  [Tig*  12). 

Again,  the  SPSS  program  FREQUENCIES  was  used  to  obtain  a 
histogram  and  the  following  descriptive  statistics:  mean, 
variance,  standard  error,  standard  deviation,  minimum  value, 
maximum  value,  and  range. 

In  Fig.  13  we  have  a  histogram  representing  the  distribu¬ 
tion  of  peak  number  flux  intensities.  The  average  value  of 

PEAK  as  determined  by  FREQUENCIES  is  9.7  x  10®  electrons/ 

2  8 
cm  -sec-ster,  with  a  standard  deviation  of  9.6  x  10 

2 

electrons/cm  -sec-ster.  The  outlyers  (enhancements  with  PEAK 

9  2 

values  greater  than  2  x  10  electrons/cm  -sec-ster)  are  the 
same  enhancements  that  are  associated  with  the  larger  mag¬ 
netic  disturbances  mentioned  in  the  previous  chapter.  Fig.  13 
also  presents  some  additional  descriptive  statistics  as  deter¬ 
mined  by  FREQUENCIES. 
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Discriminant  Analysis  of  Electron  Enhancements  and 


Their  Relationshi 

P 

to 

S 

cintillations 

The  fourth  and  last  subobjective  is  an  attempt  to  under¬ 
stand  the  relationship  between  electron  enhancement  occur¬ 
rence  and  radio  scintillations  (amplitude  and  phase).  This 
is  the  primary  objective  of  the  thesis,  and  there  are  many 
subtleties  that  one  must  understand  in  order  to  interpret  the 
results . 

The  following  procedures  must  be  used  if  one  wishes  to 
correlate  enhancements  with  scintillations.  Recall  that  the 
HILAT  vehicle  is  measuring  electrons  in  situ  at  800  km,  but 
the  scintillation  producing  irregularities  are  down  in  the  F 
region  at  350  km.  Consequently,  enhancements  must  be  spatially 
not  temporally  correlated  with  their  respective  scintillation 
values  ( S 4  and  PRMS,  with  the  two  frequencies  of  interest 
being  137  MHz  and  413  MHz). 

Using  the  17  January  pass  again  as  an  example,  and  using 
the  second  enhancement,  the  method  for  determining  the  corres¬ 
ponding  S4  and  PRMS  values  is  as  follows.  An  enlarged 
J-sensor  plot  shows  the  second  enhancement  starting  at  1:50:01 
U.T.  (Fig.  12).  One  then  tries  to  find  the  closest  corres¬ 
ponding  time  on  the  satellite  location  printout  (Table  3A) 
which  in  this  case  is  1:49:54.  This  printout  gives  the  satel¬ 
lite's  geographic  subpoint  at  this  time,  which  is  70.5  degrees 
north  latitude  and  43.9  degrees  east  longitude.  The  station 
at  Tromso  is  at  69.7  degrees  north  latitude  and  18.9  degrees 


TABLE  3-A.  Satellite  Location  Printout 


east  longitude,  so  at  this  time  the  satellite  is  east  of  the 
stat ion . 

As  we  are  investigating  the  null  hypothesis  that  enhance¬ 
ments  cause  scintillations,  the  low-energy  electrons  must 
somehow  travel  down  to  the  F  region  at  350  km.  Obviously 
they  would  have  to  spiral  about  the  magnetic  field  lines;  con¬ 
sequently,  we  need  to  find  the  latitude  where  the  field  lines 
intersect  350  km.  This  value  appears  under  the  column  headed 
FLAT  (Table  3A) ,  which  is  in  turn  under  the  heading  FLD  350  KM, 
For  this  example,  FLAT  is  at  71.3  degrees  latitude. 

We  then  look  under  the  PLAT  (penetration  latitude)  column 
of  the  scintillation  data  printout.  Table  3B,  to  find  71.3 
degrees;  71.2  (line  27)  is  what  will  be  used.  This  penetra¬ 
tion  latitude  is  also  the  latitude  at  which  the  line  of  sight 
of  the  radio  beacon  penetrates  350  km.  Note  the  two  columns 
preceding  this  value  along  line  27  are  values  of  the  topocen- 
tric  elevation  and  azimuth  of  the  350  km  penetration  point. 

For  this  example,  the  penetration  point  is  at  37.5  degrees 
elevation  and  61.0  degrees  azimuth.  (See  Fig.  14).  Note  that 
we  are  assuming  that  electron  enhancements  are  both  at  the 
same  latitude  as  the  irregularity  but  that  they  span  several 
degrees  of  longitude. 

According  to  Basu  and  Eileen  Mackenzie,  any  penetration 
point  with  elevations  less  than  20  degrees  are  suspect  and 
hence  corresponding  scintillation  values  should  not  be  used. 


Scintillation  Data  Printout 
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It  would  be  prudent  also  at  this  point  to  say  that  Dr.  Basu 
was  the  originator  of  this  procedure. 

At  this  point,  we  can  now  read  off  the  S4  and  PRMS  values 
for  the  two  frequencies  of  interest,  which  are  0.28  and  5.5 
respectively  for  137  MHz,  and  0.11  and  1.6  for  413  MHz.  The 
S4  values  are  unitless  and  the  PRMS  values  are  in  radians. 

Note  that  the  time  of  1:49:24  appears  in  the  first  column 
(U.T.)  as  the  time  of  the  previously  mentioned  scintillation 
values.  Due  to  an  error  in  formatting  the  data,  one  must  add 
7.5  seconds  to  this  time  to  get  the  actual  time  of  these  scin¬ 
tillation  values. 

On  Table  3C  we  must  perform  another  check  on  the  reliabil¬ 
ity  of  the  data  by  looking  at  the  parameter  LZEN,  or  the 
zenith  angle.  This  particular  zenith  angle  is  the  angle  be¬ 
tween  the  radio  propagation  direction  and  the  magnetic  field 
line.  Typically,  any  angle  less  than  approximately  10  degrees 
will  result  in  a  geometric  enhancement,  i.e.,  since  the 
irregularities  are  aligned  along  the  field  lines,  the  propa¬ 
gation  path  is  entering  edge  on  and  large  scintillation  can 
result  because  of  the  geometry.  Again,  for  this  example,  we 
cross-index  our  PLAT  value  (71.2  degrees)  with  the  LZEN  value, 
which  is  27.1  degrees. 

The  resulting  S4  and  PRMS  values  are  evaluated  with  re¬ 
spect  to  values  that  precede  and  follow  them,  as  we  are  inter¬ 
ested  in  seeing  if  a  rise  in  low-energy  electron  number  flux 
gives  a  corresponding  increase  in  S4  and  PRMS.  This  is 
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similar  to  the  study  of  enhancements  and  the  accompanying 
changes  in  magnetic  field.  Consequently,  we  are  interested 
in  the  rise  of  S4  and  PRMS  above  background  scintillation 
values.  Figs.  15  and  16  are  plots  of  the  S4  and  PRMS  values 
for  both  frequencies  of  this  pass,  with  the  two  scintillation 
cases  of  interest  marked  as  1  and  2.  There  are  three  enhance¬ 
ments  appearing  in  this  pass;  however,  the  LZEN  value  corres¬ 
ponding  with  the  third  one  is  10  degrees  and  hence  was  not 
used  in  this  analysis. 

For  the  purposes  of  this  research,  we  assume  that  the 
observed  enhancements  are  long-lived  (on  the  order  of  10  min¬ 
utes)  and  they  are  not  moving.  This  appears  to  be  a  reason¬ 
able  assumption,  according  to  Dr.  Hardy  (private  communica- 
t ion) . 

One  should  also  expect  that  the  observation  time  of  an 
enhancement  should  precede  that  of  the  observed  scintillation 
in  order  to  be  consistent  with  our  "cause  and  effect"  hypoth¬ 
esis.  However,  approximately  half  of  the  cases  of  interest 
exhibited  the  exact  opposite  order,  i.e.,  the  scintillation 
came  before  the  enhancement.  This  is  easily  interpreted  as 
follows.  The  satellite's  radio  beacon  has  penetrated  the  dis¬ 
torting  region  at  350  km  before  it  has  encountered  the  enhance 
ment  (the  source  of  the  irregularity)  at  the  satellite's 
orbital  altitude  of  800  km.  Further,  it  is  assumed  that  the 
electron  enhancement  has  been  going  for  some  time  (at  the 
point  in  space  of  the  satellite  encounter)  before  the  actual 
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pTfiURE  16.  Plot  of  lil-uj  for  17  January  19^4  ^ass 


satellite  encounter.  In  Fig.  17  this  situation  is  pictured 
for  the  second  enhancement  of  the  17  January  pass. 

There  should  (we  emphasize  should)  be  no  difference  be¬ 
tween  cases  where  either  the  scintillation  time  comes  before 
or  after  the  time  of  enhancement;  however,  as  a  precaution, 
a  "flag”  known  as  "T"  has  been  assigned  in  every  case.  The 
variable  T  is  either  "0"  or  "1";  if  "1",  the  scintillation  is 
observed  before  the  electron  enhancement.  If  T  is  zero,  the 
opposite  situation  exists. 

Note  also  the  difference  in  time  between  the  scintillation 
and  the  enhancement;  in  our  example  the  time  between  events  is 
23  seconds.  We  would  expect  that  for  a  small  time  difference 
we  would  best  be  able  to  determine  whether  the  enhancement's 
presence  at  800  km  affects  irregularity  development  at  350  km. 
Therefore,  another  variable  to  be  included  in  the  analysis  is 
DT,  the  difference  in  time  between  the  two  events  in  seconds. 

As  was  previously  discussed,  we  anticipate  that  for  a  rise 
in  electron  number  flux  a  commensurate  rise  in  S4  and  PRMS 
should  result.  If  a  noticeable  rise  is  evident  in  S4  and 
PRMS,  we  consider  that  a  "hit",  and  use  the  variable  "H"  where 
it  is  "1"  if  it  is  a  "hit".  A  "miss"  (no  noticeable  increase) 
has  an  H  value  of  "0".  This  particular  task  is  somewhat  easier 
said  than  done,  however,  for  the  following  reason.  Each 
scintillation  data  point  is  an  average  value  over  15  seconds; 
the  satellite  traverses  an  enhancement  in  considerably  less 
time  (on  the  order  of  three  seconds).  Consequently,  any 
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small  amplitude  or  short  time  duration  variation  in  S4  and 
PRMS  is  mostly  "buried"  in  the  averaging  process.  (The  situ¬ 
ation  is  analogous  to  photographic  resolution;  physical  fea¬ 
tures  smaller  than  the  camera's  resolution  will  not  be  easily 
detected  unless  very  bright). 

The  data  format  is  shown  in  Table  4,  where  we  have  a  total 
of  42  enhancements  with  their  associated  scintillation  values. 
Some  enhancements  here  were  not  the  same  ones  used  in  the 
previous  subobjectives  (2-4),  as  it  was  discovered  that  for 
a  few  passes  the  corresponding  scintillations  were  either 
unreliable  (low  elevation  and/or  low  LZEN)  or  defect ive .  The 
key  below  the  figure  explains  the  variable  names  used. 

The  remaining  portion  of  this  chapter  will  be  devoted  to 
explaining  discriminant  analysis  and  the  results  obtained  by 
using  this  technique.  Basically,  discriminant  analysis  is  a 
method  that  allows  one  to  statistically  distinguish  between 
two  or  more  groups  of  cases.  In  this  research,  we  will  look 
at  groups  of  scintillations  based  on  their  increase  in  S4  and 
PRMS  (three  groups;  low,  medium,  high). 

In  order  to  distinguish  between  the  groups,  a  collection 
of  discriminating  variables  that  measures  characteristics  on 
which  the  groups  may  differ  must  be  selected  (Klecka  et_  aK  , 
1975:435).  Since  we  hypothesize  that  electron  enhancements 
may  cause  scintillations,  we  have  four  variables  that  charac¬ 
terize  electron  enhancement  characteristics.  They  are  the 
peak  electron  number  flux  (PEAK),  the  time  duration  (TIME), 


average  electron  energy  (EAVG) .  These,  plus  DT,  give  a  total 


of  five  discriminating  variables. 

The  discriminating  aspect  is  accomplished  by  weighting 
and  then  linearly  combining  the  discriminating  variables  in 
order  to  make  the  groups  as  distinct  as  possible.  These  "dis¬ 
criminant  functions"  are  of  the  form 


D i  ~  dilZl  +  d i 2 Z 2  + 


+  d.  Z 

ip  p 


where  Di  is  the  score  of  the  discriminant  function  i,  the  d's 
are  the  weighting  coefficients,  and  the  Z's  are  the  standard¬ 
ized  values  of  the  p  discriminating  variables.  The  maximum 
number  of  derivable  functions  is  either  one  less  than  the  num¬ 


ber  of  groups  or  equal  to  the  number  of  discriminating  vari- 


particular  discriminating  variable  contributes  to  the  dif¬ 
ferentiation  (Klecka  et  al . ,  1975:436). 


The  classification  aspect  arises  after  the  discriminating 
functions  have  been  computed.  After  the  necessary  discrimina¬ 
tion  variables  have  been  determined  for  cases  with  known  group 
membership,  classification  functions  can  be  calculated  which 
will  allow  future  classification  of  cases  of  unknown  member¬ 
ship.  The  equation  for  a  particular  group  would  appear  as 


‘  cilVl  *  ci2V2  * 
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where  C ^  is  the  classification  score  for  group  i,  the  c's  are 
the  classification  coefficients,  and  the  V's  are  the  variable 
values  of  the  discriminating  variables  (Klecka  et.  jd. ,  1975: 
436,445) . 

Under  the  assumption  of  a  multivariate  normal  distribution 
classification  scores  can  be  used  to  determine  probabilities 
of  group  membership.  In  other  words,  by  assigning  a  case  to 
the  group  with  the  highest  value  of  V,  we  are  saying  that 
this  case  has  the  greatest  probability  of  being  in  group  i 
(Klecka  et  al.  ,  1975:445). 

The  value  of  the  classification  functions  are  to  tell  us 
how  effective  the  discriminating  variables  are.  If  there  is 
a  great  number  of  misclass  if icat ions ,  then  more  reliable  dis¬ 
criminating  variables  need  to  be  found.  The  routine  used  in 
this  research  allows  for  a  display  for  each  case  of  the  dis¬ 
criminating  score  and  the  predicted  group  membership.  The 


largest  probability  is  given,  followed  by  the  next  highest 
probability  if  larger  than  0.0005.  Another  probability  given 
is  the  probability  that  a  member  of  the  predicted  group  would 
be  as  far  from  the  group  centroid  (the  mean  value  of  a  dis¬ 
criminating  variable  averaged  over  all  groups)  as  the  case 
under  consideration;  this  probability  is  represented  as  P(X/G) 
(Klecka  et  al.  ,  1975:446). 

Two  methods  were  used  to  determine  the  desired  discrimina¬ 
tion.  The  first  method  is  direct,  i.e.,  all  discriminating 
variables  are  used.  In  the  second  method,  a  stepwise  selec¬ 
tion  procedure  is  used.  Based  on  a  user-defined  criterion, 
the  procedure  starts  with  the  single  best  discriminating  vari¬ 
able.  A  second  discriminating  variable  is  selected  if  when 
combined  with  the  first  variable  greater  discrimination  results 
Third  and  subsequent  variables  are  selected  at  each  remaining 
step,  along  with  simultaneous  variable  elimination  if  greater 
discrimination  results.  After  all  variables  have  been  con¬ 
sidered  (or  if  the  remaining  variables  do  not  render  increased 
discrimination),  the  procedure  stops  (Klecka,  1975:436). 

In  this  case,  the  user-defined  criterion  was  to  maximize 
the  smallest  Mahalonobis  distance  (McNichols,  1980:(7-44)- 
(7-48))  between  two  groups.  The  Mahalonobis  distance  can  best 
be  defined  by  starting  with  the  definition  of  Euclidean  dis¬ 
tance  in  p-d  imens ional  space.  Two  points  x.  and  x9  are  a 

X  u 

2 

squared  distance  D  apart;  in  vector  terms  this  can  be  ex¬ 


pressed  as  (Fig.  18) 


^TBIJRE  IB.  Euclidean  Distance  ( McNi ehols, 1 9B0: ( 7-44 ) ) 


D  =  (x1-x2)'(x1-x2)  =  EU1.-x2j)  (28) 

Intuitively,  one  would  expect  a  case  that  is  closest  to  a 
particular  group  centroid  (smallest  D)  would  be  assigned  to 
that  group.  The  problem  with  this  approach  can  be  seen  by 
examining  Fig.  19  that  represents  a  multivariate  population 
distribution,  looking  down  on  an  x^-x2  plane  in  the  case  with 
two  variables  (McNichols,  1980: (7-45) ) . 

The  reason  why  Euclidean  distance  is  not  used  has  to  do 
with  the  fact  that  the  distance  along  the  x2  dimension  is  more 
meaningful  statistically  than  in  the  x^  dimension;  this  is 
because  the  x^  standard  deviation  is  larger  than  the  x2  stan¬ 
dard  deviation.  Therefore,  an  observation  has  a  greater 
possibility  of  being  one  unit  away  from  the  centroid  in  the 
x1  direction  than  in  the  x2  direction  (McNichols,  1980: (7-45)- 
(7-46)). 

In  order  to  compensate  for  this,  both  the  x^  and  x2 
coordinates  are  divided  by  their  respective  standard  devia¬ 
tions.  With  this  correction,  the  squared  distance  between 
any  observation  x^  and  a  centroid  x^  will  be 

1)2  =j,[(xu-'xkj)/0jl2  (29) 

Rewritten  in  matrix  notation 
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expression  as  £  ^ ,  we  finally  have 

D2  =  (xi-xk)'£"1  (xrxk)  (31) 

2 

This  final  expression  for  D  is  then  the  square  of  the 
Mahalonobis  distance  (McNichols,  1980 : (7- 46) - (7-48) ) . 

How,  then,  does  one  evaluate  the  printed  output  generated 
by  the  discriminant  analysis  subprogram?  The  subprogram  used 
(DISCRIMINANT)  outputs  analytic  features  that  allow  one  to 
evaluate  the  importance  of  the  computed  discriminant  functions, 
as  well  as  how  to  interpret  the  computed  discriminant  function 
coefficients.  Discussion  of  these  features  appear  in  the 
ensuing  paragraphs. 

Space  does  not  allow  a  full  mathematical  development  of 
discriminant  analysis,  but  the  two  main  assumptions  that  must 
be  presented  are  that  the  discriminating  variables  have  a 
multivariate  normal  distribution,  and  that  they  have  equal 
variance-covariance  matrices  within  each  group.  Fortunately, 
the  technique  is  sufficiently  robust  to  allow  some  departure 
from  the  stated  assumptions  (Klecka,  1975:435).  At  any  rate, 
in  order  to  determine  the  coefficients,  an  eigenvector  prob¬ 
lem  of  the  type 


WA  =  A  Ba  (32) 

must  be  solved.  Here  B  and  W  are  the  "between-groups"  and 
"within-groups"  sums  of  squares  and  crossproducts,  respec¬ 
tively. 
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For  each  discriminant  function,  there  is  a  corresponding 
eigenvalue  X  .  Since  there  may  be  several  functions,  it  would 
be  desirable  to  assess  the  importance  of  one  over  the  others. 
The  eigenvalue  is  such  a  measure;  it  is  the  relative  impor¬ 
tance  of  each  function,  where  the  sum  of  the  eigenvalues  is 
a  measure  of  the  total  variance  existing  in  the  discriminating 
variables  (Klecka,  1975:442). 

Another  aid  in  judging  the  importance  of  a  discriminating 
function  is  its  canonical  correlation.  The  canonical  correla¬ 
tion,  when  squared,  is  the  proportion  of  variance  in  the  dis¬ 
crimination  function  as  explained  by  the  groups  (Klecka,  1975: 
442)  . 

A  second  method  used  to  determine  importance  is  to  test 
for  the  statistical  significance  of  discriminating  information 
not  already  incorporated  in  earlier  functions.  Starting  with 
no  functions,  and  for  each  function  thereafter,  Wilk's  lambda 
is  computed.  Wilk's  lambda  is  an  inverse  measure  of  the  dis¬ 
criminating  power  left  in  the  original  variables  not  yet  re¬ 
moved  by  the  discriminating  functions.  It  can  also  be  convert¬ 
ed  into  a  chi-square  statistic  in  order  to  test  for  statisti¬ 
cal  significance  (Klecka,  1975:442-443). 

As  was  mentioned  previously,  the  coefficients  are  a 
measure  of  each  discriminating  variable's  contribution  to  the 
resulting  discriminating  function.  DISCRIMINANT  yields  both 
standardized  and  unstandardized  coefficients;  standardized 
coefficients  are  derived  so  that  the  resulting  discriminant 
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score,  Di>  is  in  standard  form.  This  means  that  for  all  cases 
used  in  the  analysis,  the  score  from  one  function  will  have  a 
mean  of  zero  and  a  standard  deviation  of  one  (Klecka,  1975: 

443) . 

For  computational  purposes,  standardized  coefficients  would 
be  of  little  use  as  the  discriminating  variables  are  rarely 
in  standard  form;  therefore,  unstandardized  coefficients  are 
given  if  they  are  desired.  However,  unstandardized  coeffi¬ 
cients  cannot  be  used  to  assess  relative  importance  because 
they  lack  adjustment  for  both  measurement  scales  and  the 
variability  found  in  the  original  variables  (Klecka,  1975: 
443-444)  . 

After  performing  some  initial  runs,  it  became  evident  that 
some  consistent  method  was  needed  to  classify  the  rise  in 
scintillation  values  (DS1,  DS2,  DPI,  DP2)  in  order  to  yield 
reliable  results.  Here  DS1  is  the  rise  in  S4  at  137  MHz,  and 
DS2  is  the  rise  in  S4  at  413  MHz.  Similarly,  DPI  is  the  rise 
in  PRMS  at  137  MHz  and  DP2  is  the  rise  in  PRMS  at  413  MHz. 

By  running  FREQUENCIES,  which  prints  out  a  histogram,  it  was 
easier  to  get  an  idea  of  how  the  low,  medium,  and  high  values 
of  the  rise  in  scintillation  were  distributed.  A  histogram 
for  DS1  appears  in  Fig.  20;  DS1,  DS2,  DPI,  and  DP2  were  broken 
out  in  terms  of  low,  medium,  and  high  values.  Here  the  pri¬ 
mary  criteria  for  division  was  to  ensure  as  near  an  equal 
size  as  possible  of  the  three  groups.  Also,  graphs  of  DS1 
versus  DS2  (Fig.  21)  and  DPI  versus  DP2  (Fig.  22)  are  shown. 
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[G-URE  20-D.  Histogram  of  DP2  (42  cases) 


It  was  quite  difficult  to  classify  the  DS1  and  DS2  values 
because  they  often  did  not  increase  or  decrease  in  concert 
with  each  other;  the  DPI  versus  DP2  graph  on  the  other  hand 
gives  a  more  consistent  pattern,  i.e.,  for  a  rise  in  DPI  there 
is  a  commensurate  rise  in  DP2. 

This  fact  is  significant  as  it  validates  the  expression 
(Eq  (1))  for  PRMS  as  put  forth  by  Briggs  and  Parkin.  Using 
this  equation,  we  can  establish  the  following  ratio  for  DPI 
and  DP2 

DP1/DP2  =  > ,  /  X  ^  (33) 

as  the  other  variables  in  Eq  (1)  are  the  same  for  the  two 
radio  frequencies.  In  this  analysis,  X,  =  2.2  m  (since  f1  =  137 
MHz)  and  Aa  =  .73  m  (since  f2=413  MHz),  making  X,  /  X  ^ 
approximately  3.  In  Fig.  23,  DP2  is  on  the  vertical  axis  and 
DPI  is  on  the  horizontal  axis,  so  that  the  slope  of  the  line 
is  about  .33,  or  Xa/X,as  shown.  Notice  that  the  data  agrees 
fairly  well  with  the  theory. 

Briggs  and  Parkin  presented  expressions  for  S4  that  did 
not  show  such  a  linear  dependence  on  as  in  the  case  of  PRMS; 
this  may  be  the  reason  for  the  more  scattered  appearance  of 
the  distribution  for  DS1  versus  DS2  as  exhibited  in  Fig.  21. 

We  might  expect  that  the  data  would  fall  between  the  two  curves 
that  appear  in  Fig.  21;  these  curves  represent  the  limiting 
cases  of  Eq  (10).  The  upper  line  is  Eq  (11), 


S(X,)/sog  =  X,/X> 


(ID 


AD-A162  839  STATISTICAL  ANALYSIS  OF  LOU  ENERGY  ELECTRONS  AND  THEIR  2/2 
RELATIONSHIP  TO  RA.  .  <U)  AIR  FORCE  INST  OF  TECH 
HRIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENGI. .  M  E  OBER 
UNCLASSIFIED  DEC  84  AFIT/GS0/ENP-ENS/B4D-2  F/G  4/1  NL 


as  applied  to  this  particular  set  of  data  (small  value  of 

r  ,  the  traverse  radius  of  the  irregularity),  and  the  lower 
o 

line  represents  the  case  where  rQ  is  large  (Eq  (12)). 

S(\)/S(X>)  =  (*,/>a)2  (12) 

We  must  also  analyze  Eq  (10)  in  order  to  determine  whether 
or  not  it  is  a  monotonically  increasing  or  decreasing  func¬ 
tion  of  r  .  This  is  a  necessary  condition  if  all  S4  data 
o 

points  are  to  be  bounded  by  the  curves  denoting  large  and  small 

values  of  r  .  This  can  be  done  by  differentiating  Eq  (10) 
o 

with  respect  to  r  ,  and  then  set  the  derivative  to  equal  zero. 
If  there  is  a  value  of  r0 that  makes  Eq  (10)  a  minimum  or  maxi¬ 
mum,  then  we  can  say  that  it  is  not  a  monotomic  function. 

Eq  (10)  is  repeated  here  for  easier  reference, 

-  */x  ‘/a. 

S(\)/S(XJ)  =A,/AJ.(l  +T*V/HXla23)(l+T»V/^^2?)  (10) 

Differentiating  and  setting  equal  to  zero  using  the  product 
rule  yields 


x.A.C-i  (I  t-oV.*  (V«i  a‘zh(  I  ftrV  /H  x*zY+ 

(  V*  CirV/X‘£ L)0-tirO /t\^z  )  =  O 


Simplification  leave 


The  only  way  that  Eq  (35)  can  be  true  is  if  rQ  is  zero 
id  clearly  A, is  not  equal  to  A^in  this  case.  As  a 


* .  •  _ 
r.V.  -\ 

J|V  ,1*l1 
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further  check,  Eq  (10)  is  plotted  for  a  range  of  rQ  values 
(1-4000  m,  Fig.  23).  The  value  Z,  which  is  the  "reduced  dis¬ 
tance"  1^2!  {z\tZ2 )  comPuted  for  the  case  when  the  satel¬ 
lite  is  directly  overhead.  Taking  the  height  of  the  irregu¬ 
larity  =  350  km,  %2  is  then  the  satellite  altitude  h  = 

800  km  minus  z^.  This  gives  a  Z  of  approximately  200  km. 

In  Fig.  23  and  in  Eq  (10),  the  limiting  case  S(>,)/S(A)  =  \l\x 
results  when  TJ  rcH/4  N*  Z*<<  1  andir1r'/4>fz*‘ <<  1.  The  other 
limiting  case  S(>,)/S(Xl)  =  (X,/^)^  results  when  if*  reH/  4  X*  Z*-  >  >  1 
and  If  r*/4X*Z  >>  1.  In  the  first  case,  if  Tf3  rj4  /4  X*  Z \  0 . 1 , 
then  rQ  is  on  the  order  of  10  m.  In  the  second  case,  if 
If  r„/4X  Z  >10,  then  rQ  is  on  the  order  of  a  kilometer. 

A  significant  amount  of  the  data  appearing  in  Fig.  21  lie 
in  the  vicinity  of  the  two  curves.  There  are  a  few  outliers 
(DS1  less  than  .10,  DS2  less  than  .05);  however,  the  reason 
for  this  is  not  clear. 

One  point  of  concern  dealt  with  the  independence  of  the 
discriminating  variables  with  respect  to  each  other;  in  par¬ 
ticular,  the  time  integrated  number  flux  PTIME.  The  reason 
being  that  for  each  case  PTIME  was  computed  from  PEAK  and 
TIME.  For  example,  in  the  case  where  the  area  under  the  PEAK 
curve  approximated  a  square  PTIME  was  just  PEAK  times  TIME. 

In  the  case  where  the  area  under  the  PEAK  curve  was  not  so 
clear  cut,  it  was  subdivided  into  simple  square  and  triangular 
areas  so  that  PTIME  was  basically  the  summed  area  of  the 
smaller  subdivisions,  still  making  it  a  function  of  PEAK  and 
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Fortunately,  a  Pearson  correlation  coefficient  routine  is 
available  in  SPSS;  the  Pearson  correlation  coefficient  is  a 
measure  of  the  strength  and  direction  (negative  being  inversely 
proportional,  positive  being  directly  proportional)  of  corre¬ 
lation  between  two  variables.  Table  5  contains  the  results 
of  pairwise  comparisons  of  the  five  discriminating  variables; 
note  that  both  PEAK  and  TIME  are  highly  correlated  with  PTIME 
(.7087  and  .7064,  respectively).  Two  series  of  discriminant 
analysis  runs  were  subsequently  made;  one  with  PTIME  and  one 
without  PTIME  for  the  sake  of  comparison. 

In  Table  6  we  have  the  results  of  a  stepwise  selection 
discriminant  analysis  run,  where  the  variable  CAT  is  the  over¬ 
all,  or  average  score  for  a  particular  case  based  on  all  four 
scintillation  variables.  This  is  done  by  using  the  classifi¬ 
cation  histograms  (Fig.  20);  for  a  given  case,  the  classifi¬ 
cation  for  each  variable  (DS1,  DS2,  DPI,  and  DP2)  is  noted. 
Letting  low=l,  medium=2,  and  high=3,  the  overall  classifica¬ 
tion  is  obtained  by  adding  up  the  four  values  and  then  divid¬ 
ing  by  four.  The  problem  with  this  approach  is  that  for  a 
given  case  one  might  have  a  low  DS1  value,  a  medium  DS2  value, 
a  high  DPI  value,  and  a  medium  DP2  value;  it  became  evident 
that  this  method  failed  to  give  meaningful  results. 

Instead  of  attempting  a  discriminant  analysis  on  the  42 


cases  of  interest  based  on  overall  classification  for  each 
case,  eight  different  types  of  classifications  were  devised. 
They  are  C137  (DS1  and  DPI),  C413  (DS2  and  DP2 ) ,  CS4  (DS1 


TABLE  5.  Pearson  Correlation  Coefficient  Printout 
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Table  6.  Classification  Based  on  CAT,  w/PTIME 


After  Function 

Wilk ' s \  Sig 

Funct ion 

Eigenvalue 

0 

.766  .119 

1 

.  168 

1 

.895  .121 

2 

.118 

Standardized 

Discr iminat ion 

Coefficient  s 

Fund 

Func2 

TIME 

.44469 

1.35076 

PTIME 

.24721 

-1.25659 

DT 

- .65380 

.27246 

Classification  Results 

Actual  Group 

#  of  cases 

Predicted  Group 
1  2 

Membership 

3 

Group  1 

13 

8  3 

(61.5)  (23.1) 

2 

(15.4) 

Group  2 

15 

5  7 

(33.3)  (46.7) 

3 

(20.0) 

Group  3 

14 

4  3 

(28.6)  (21.4) 

7 

(50.0) 

%  of  grouped  cases  correctly  classified  -  52.38 


and  DS2) ,  CPRMS  (DPI  and  DP2) ,  CS1  (DS1) ,  CS2  (DS2) ,  CPI 
(DPI),  and  CP2  (DP2).  Again,  for  each  of  the  eight  types 
of  classification,  the  same  three  groups  based  on  the  histo¬ 
grams  appearing  in  Fig.  20  are  used.  For  C137,  C413,  CS4, 
and  CPRMS,  overall  classification  is  decided  by  summing  the 
value  of  each  variable  (again  from  Fig.  20)  and  then  dividing 
by  two;  in  the  event  that  a  non-integer  result  was  obtained, 
the  overall  classification  was  determined  by  rounding  down. 

Another  concern  that  arose  out  of  the  initial  discriminant 
analysis  runs  was  that  many  of  the  observed  low  changes  in 
scintillation  occurred  before  an  enhancement  was  encountered, 
i.e.,  T=l.  When  no  scintillation  occurs,  it  is  conceivable 
that  the  observed  enhancement  may  have  begun  just  at  the  time 
it  is  observed  by  the  satellite.  Hence,  a  series  of  eight 
runs  (one  for  each  type  of  classification)  without  the  ques¬ 
tionable  cases  (T=l  and  low  scintillation  value)  was  per¬ 
formed.  The  total  number  of  cases  then  were  as  high  as  38, 
with  a  low  of  30.  For  eight  runs,  the  average  number  of 
cases  was  approximately  33.  Low  scintillation  values  for  CS1 
and  CS2  were  DS1  and  DS2  values  less  than  or  equal  to  .02. 

Low  values  for  CPI  were  DPI  values  less  than  or  equal  to  .4, 
and  low  CP2  values  were  less  than  or  equal  to  .2  for  DP2. 

The  low  values  for  C137  (DS1*.02  and  DPI*. 4),  C413  (DS2*.02 
and  DP24.2),  CS4  (both  DS1  and  DS24.02),  and  CPRMS  (DPI*. 4 
and  DP2*.2)  are  based  on  the  same  criterion  for  low  scintilla¬ 
tion  values  as  CS1,  CS2,  CPI,  and  CP2. 


In  Table  7  we  have  the  results  of  the  discriminant  analy¬ 
sis  based  on  CPRMS  without  the  use  of  PTIME.  Two  functions 
were  computed,  with  the  first  being  more  important  than  the 
second  based  on  its  larger  eigenvalue  (.361  compared  to  .091). 
In  the  upper  right  hand  corner,  notice  that  a  good  deal  of 
discriminating  power  exists  in  the  first  function  as  evidenced 
by  the  Wilk's  lambda  value  of  .674.  In  this  research  effort, 
a  Wilk's  lambda  with  a  corresponding  significance  value  less 
than  .10  is  considered  to  be  statistically  meaningful.  In 
this  example,  the  computed  significance  is  .020,  so  that  the 
first  discriminant  function  will  yield  separation  between  the 
groups  that  will  not  be  solely  due  to  chance.  The  second 
function  has  a  corresponding  Wilk's  lambda  of  .916  and  a  sig¬ 
nificance  of  .191,  indicating  that  it  will  not  aid  much  in 
providing  discrimination  between  the  groups. 

Upon  inspection  of  the  computed  discriminant  function 
coefficients  for  both  functions,  we  see  that  in  Function  1 
TIME  and  DT  make  the  largest  contribution  (sign  indicating 
positive  or  negative  contribution),  with  EAVG  being  of  less 
importance.  In  Function  2,  which  is  not  significant,  EAVG 
makes  the  largest  contribution  followed  by  DT,  with  TIME  being 
the  lesser  contributor  of  the  three. 

Based  on  the  computed  discrimination  functions,  the 
classification  results  appear  at  the  bottom  of  the  table. 

For  example,  there  were  19  Group  1  members  (low  scintillation) 
to  start  with;  the  computed  classification  functions  placed 


Table  7.  Classification  Based  on  CPRMS  w/o  PTIME 


After  Function 


Funct ion 

Eigenvalue 

0 

1 

.361 

1 

2 

.091 

Standardized 

Discr iminat ion 

Coef f ic ients 

Func  1 

Func2 

TIME 

.82785 

.19298 

EAVG 

-.02290 

.95978 

DT 

- .46279 

.55163 

Classification  Results 

Actual  Group 

#  of  cases 

Predicted  Group 
1  2 

Group  1 

19 

13  5 

(68.4)  (26.3) 

Group  2 

12 

3  8 

(25.0)  (66.7) 

Group  3 

11 

4  4 

(36.4)  (36.4) 

%  of  grouped  cases  correctly  classified  - 


Wilk ' sX  Sig 
.674  .020 
.916  .191 


Membership 

3 

1 

(  8.3) 

1 

(  8.3) 

3 

(27.3) 

57.14 


13  of  the  Group  1  members  in  Group  1,  five  in  Group  2,  and 
one  in  Group  3.  The  Group  2  members  were  not  classified  quite 
as  well,  and  the  Group  3  members  were  for  the  most  part  in¬ 
correctly  classified.  The  overall  percentage  of  cases  cor¬ 
rectly  classified  is  57.14,  or  about  24  of  the  42  cases. 

Another  similar  discriminant  analysis  based  on  PRMS  was  made 
that  included  PTIME  as  a  discriminating  variable;  however, 
the  final  results  were  identical  to  those  presented  above. 

The  other  three  classifications  that  rendered  signifi¬ 
cant  results  were  based  on  C413,  CS1,  and  CP2  (see  Appendix  C 
and  Table  8) .  Note  that  in  the  CS1  run  there  was  a  difference 
in  final  results  when  PTIME  was  included;  however,  PTIME' s 
inclusion  does  little  to  increase  separation  between  the 
groups.  A  similar  situation  exists  with  C413;  without  PiIME 
the  results  are  significant,  but  with  PTIME's  inclusion  both 
of  the  eigenvalues  and,  therefore,  the  discriminant  functions 
are  not  significant.  Table  8  contains  a  summary  of  the 
results  of  all  eight  runs  where  all  42  cases  are  used. 

In  the  other  16  runs  where  the  questionable  low  valued 
observations  were  discarded,  the  classifications  that  rendered 
significant  results  were  C413  and  CPRMS  (see  Appendix  D) . 

The  overall  results  for  this  series  of  runs  appears  in  Table  9. 

In  Table  10  we  have  a  comparison  of  the  two  series  of 
runs.  For  each  classification,  the  significance  for  each 
series  (all  42  cases  versus  the  situation  where  the  question¬ 
able  cases  are  eliminated,  designated  T=l)  and  the 
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Table  8. 

Discriminant  Analysis 
(both  T=0  and 

Results 
T=  1) 

For  All  Cases 

Number  of 

After 

Cases 

Classification  Function 

Wilk ' 

'  s  X  Sig 

42 

C137 

0 

1 

.929 

.240 

42 

*C413(w/oPTIME) 

0 

1 

.858 

.051 

42 

C413 (w/PTIME) 

0 

.830 

.  127 

1 

.996 

.678 

42 

CS4 

0 

1 

.910 

.159 

42 

*CPRMS 

0 

.674 

.020 

1 

.916 

.191 

42 

*CS1 (w/oPTIME) 

0 

.  798 

.069 

1 

.916 

.067 

42 

*CS1 (w/PTIME) 

0 

.  749 

.089 

1 

.872 

.074 

42 

CS2 

0 

(no  ! 

results) 

1 

42 

CPI (w/oPTIME) 

0 

1 

.945 

.329 

42 

CPI (w/PTIME) 

0 

.843 

.  161 

1 

.966 

.251 

42 

*CP2 

0 

.  700 

.035 

1 

.921 

.  208 

*  S i gn if  icant 
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Table  9.  Discriminant  Analysis  Results  For 
Only  T=1  Cases 


Number  of 
Cases 


Classif icat ion 


Alter 
Funct ion 


Wilk'sX  Sig 


1 

34 

C137  i 

0 

1 

.921  .277 

32 

*C4 13  i 

0 

1 

.679  .026 

.883  .059 

1 

30 

CS4 

0 

1 

.901  .244 

33 

*CPRMS 

0 

1 

> 

34 

CS1 

0 

1 

.764  .234 

.903  .217 

31 

CS2 

9 

(no  results) 

i 

38 

CPI (w/oPTIME) 

0 

1 

38 

CPI (w/PTIME) 

.750  .135 

.945  .383 

1 

34 

CP2 

.871  .118 

Signif icant 


H] 


discriminating  variables  that  appear  in  the  discriminant 
function  are  shown.  Next  to  the  discriminating  variables 
appears  the  corresponding  discrimination  coefficients  rounded 
off  to  two  significant  digits  for  comparison  purposes. 

The  results  can  be  interpreted  in  light  of  the  results 
presented  by  Briggs  and  Parkin.  The  classification  based  on 
CS1  (rise  in  S4  at  137  MHz)  shows  that  PEAK  (electron  number 
flux  peak  value)  and  EAVG  (average  energy)  are  important  dis¬ 
criminators.  Consider  the  situation  where  the  significance 
level  is  .067  (PTIME  not  included).  Upon  inspection  of  Eq  (5) 
we  see  that  S4  is  proportional  to  PRMS  which  is  in  turn  pro¬ 
portional  to  electron  number  density  N  (Eq  (1)).  This  is  con¬ 
sistent  with  the  fact  that  PEAK  is  a  good  discriminating  vari¬ 
able.  EAVG  has  a  large  positive  contribution  as  well  (.64), 
indicating  that  the  higher  the  energy  of  the  electron,  the 
greater  the  collisional  ionization  in  the  ionosphere.  This  is 
again  a  condition  one  would  expect  necessary  to  sustain  irreg¬ 
ularity  production,  which  in  turn  would  increase  S4  scintilla¬ 
tion.  Since  amplitude  scintillations  at  137  MHz  are  larger 
than  S4  values  at  413  MHz,  such  a  dependence  on  PEAK  and  EAVG 
would  be  more  evident  for  the  CS1  case;  this  may  explain  why 
a  similar  result  did  not  appear  in  the  CS2  runs. 

One  would  expect  similar  results  in  the  case  of  CPI  and 
CP2  (category  based  on  the  rise  in  PRMS  at  137  MHz  and  413 
MHz,  respectively),  yet  only  for  the  CP2  classification  do  we 
get  significant  results.  Since  PRMS  is  linear  in  number 


density  N,  we  might  expect  the  same  discriminant  variables  as 
in  the  case  of  CS1.  However,  the  discriminant  variables  that 
gave  the  best  discrimination  for  classification  based  on  CP2 
were  TIME  (-.65),  DT  (.70),  and  EAVG  (.14). 

In  the  case  where  classification  is  based  on  both  DPI  and 
DP2  together  (CPRMS) ,  we  get  similar  results:  TIME  (.83), 

DT  (-.46),  and  EAVG  (-.02).  In  both  cases  (CP2  and  CPRMS), 
the  physical  size  of  the  enhancement  and  the  time  between 
enhancement  encounter  and  scintillation  occurrence  are  impor¬ 
tant  indicators.  In  each  instance,  the  average  energy  EAVG 
plays  a  minor  role  in  discrimination. 

The  only  other  case  that  gave  statistically  meaningful 
results  was  C413;  since  its  discriminating  variables  are  both 
TIME  and,  in  the  case  where  questionable  low  scintillations 
were  discarded,  DT  (.99).  This  result  is  similar  to  that 
found  with  the  CPRMS  and  CP2  classifications.  Since  the  CS2 
classification  failed  to  produce  any  results  whatsoever,  we 
conclude  that  the  significant  discrimination  in  the  case  of 
C413  is  most  likely  due  to  the  significant  discrimination  of 
DP2.  With  this  information,  the  discriminant  analysis  based 
on  C137  and  C413  are  probably  redundant. 


VIII.  Conclusions  and  Recommendations 

Conclusions 

The  198  low  energy  electron  enhancements  were  found  to 
exist  primarily  in  the  nighttime  sectors,  with  a  large  number 
occurring  in  both  the  0300  sector  and  in  the  2300  sector. 

There  was  a  sharp  dropoff  in  enhancements  shortly  after  0700 
MLT,  with  none  (except  for  one  instance  at  0800  MLT)  appear¬ 
ing  until  approximately  1500  MLT. 

Of  the  73  enhancements  that  were  used  to  study  the  possi¬ 
bility  of  associated  Birkeland  currents,  only  nine  enhance¬ 
ments  had  such  corresponding  currents.  It  appeared  that  there 
may  be  a  correlation  between  the  number  flux  (PEAK)  and  the 
computed  current  density  J„ ;  more  statistics  are  needed  for 
such  a  study. 

O 

The  average  value  of  PEAK  was  found  to  be  9.7  x  10 
2 

electrons/cm  -sec-ster,  with  a  wide  standard  deviation  of 
8  2 

9.6  x  10  electrons/cm  -sec-ster.  Enhancements  with  PEAK 

9  2 

values  greater  than  2  x  10  electrons/cm  -sec-ster  were  the 
same  ones  that  had  corresponding  geomagnetic  field  distur¬ 
bances  as  mentioned  above. 

The  plot  of  the  PRMS  data  (DPI  versus  DP2)  agreed  well 
with  the  theory  as  presented  by  Briggs  and  Parkin,  i.e.,  PRMS 
is  linear  in  A  .  The  plot  of  the  S4  data  (DS1  versus  DS2)  was 
not  nearly  as  linear;  after  plotting  curves  for  the  limiting 
cases  of  the  equation  for  the  ratio  of  DS1  to  DS2,  many  data 
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points  did  fall  in  the  vicinity  of  the  two  curves;  however, 
many  did  not.  Why  this  is  so  is  not  clear. 

After  performing  several  discriminant  analysis  runs,  the 
/ 

results  that  were  statistically  meaningful  (for  42  cases) 
appeared  in  the  runs  where  classification  was  based  on  C413, 
CPRMS,  CS1,  and  CP2.  When  we  eliminated  low  scintillation 
cases  where  the  scintillation  occurred  before  the  enhancement 
encounter  (T=l) ,  the  classifications  that  gave  significant 
results  were  based  on  C413  and  CPRMS, 

In  the  situation  where  the  discriminant  analysis  was  based 
on  CS1  (classified  by  rise  in  S4  at  137  MHz),  the  discrimi¬ 
nating  variables  were  PEAK  and  EAVG.  Since  S4  is  proportional 
to  PRMS,  which  is  in  turn  proportional  to  number  density  N, 
then  the  fact  that  PEAK  is  a  discriminating  variable  is  con¬ 
sistent  with  theory.  We  also  expect  EAVG  to  be  an  indicator 
of  amplitude  scintillation  as  well;  higher  electron  energies 
will  result  in  increased  collisional  ionization.  This,  in 
turn,  will  contribute  to  irregularity  development,  which  is 
responsible  for  amplitude  scintillation.  Since  DS1  (rise  in 
S4  at  137  MHz)  is  larger  than  DS2  (rise  in  S4  at  413  MHz), 
this  may  explain  why  discrimination  analysis  was  successful 
in  the  former  case  and  not  the  latter. 

In  the  discriminant  analysis  results  where  the  size  of 
PRMS  was  the  criteria  of  classification,  CP2  (classification 
based  on  the  rise  in  PRMS  at  413  MHz)  and  CPRMS  (classifica¬ 
tion  based  on  the  rise  in  both  PRMS  values)  both  had  TIME 
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(time  necessary  for  the  satellite  to  traverse  the  enhancement, 
i.e.,  its  physical  size)  and  DT  (time  difference  between 
scintillation  occurrence  and  enhancement  encounter;  this  is 
indicative  of  the  length  of  the  radio  beam  propagation  path) 
as  the  major  discriminating  variables.  In  both  situations, 
EAVG  played  only  a  minor  role  in  discrimination. 

The  fact  that  classification  by  C413  (S4  and  PRMS  at  413 
MHz)  produced  significant  results  is  probably  due  to  DP2's 
ability  to  be  grouped  significantly.  Since  the  discriminat¬ 
ing  variables  are  the  same  for  this  situation  as  for  the  CP2 
case,  the  discriminant  analysis  based  on  C413  (and  C137)  is 
probably  redundant. 

Recommendat ions 

As  was  mentioned  in  Chapter  VII,  the  scintillation  data 
points  used  in  this  thesis  effort  were  averaged  over  15  sec¬ 
onds,  thereby  making  it  far  more  difficult  to  establish  a 
relationship  between  the  enhancements  and  observed  scintilla¬ 
tions,  especially  if  the  latter  are  small.  A  way  to  allev¬ 
iate  this  situation  would  be  to  increase  the  time  resolution 
of  the  scintillation  data  to  something  on  the  order  of  one 
second  as  opposed  to  fifteen  seconds.  Whether  this  can  easily 
Le  accomplished  with  the  HILAT  satellite  is  unknown;  any 
future  spaceborne  radio  scintillation  experiment  attempting 
to  investigate  the  enhancement /sc int i 1 lat ion  relationship 
might  consider  this  recommendation. 
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A  second  recommendation  is  that  at  least  an  order  of 
magnitude  more  enhancements  need  to  be  analyzed  in  order  to 
improve  statistical  results.  Time  did  not  allow  all  three 
months  worth  of  data  to  be  analyzed  in  any  detail;  all  passes 
during  the  three  month  period  that  contain  enhancements  should 
be  studied.  In  fact,  a  study  done  with  six  months  worth  or 
more  of  HILAT  data  might  give  more  conclusive  results. 

This  would  allow  one  to  gather  more  data  on  magnetic  field 
disturbances,  and  by  performing  a  linear  regression  analysis, 
more  could  be  learned  about  the  relationship  between  enhance¬ 
ments  and  magnetic  disturbances.  It  would  also  be  interesting 
to  see  if  enhancements  correlate  with  both  magnetic  field  dis¬ 
turbances  and  scintillations;  this  was  an  early  considera¬ 
tion  in  this  thesis,  but  a  lack  of  data  prevented  such  an  anal¬ 
ysis.  More  enhancement  data  would  also  allow  for  more  accu¬ 
rate  descriptors  such  as  the  average  peak  number  flux. 

The  discriminant  variables  describing  the  enhancements  in 
this  analysis  may  not  be  the  best  ones  available  to  explain 
how  scintillations  change.  There  may  indeed  be  some  charac¬ 
teristic  inherent  in  enhancements  that  render  better  separa¬ 
tion  than  those  used  here. 

A  final  recommendation  is  a  long  time-period  study  to 
track  the  effect  of  solar  activity  on  scintillation  occurrence. 
Although  traditionally  Kp  has  been  used  as  an  indicator  of 
solar  activity,  a  better  indicator  is  solar  wind. 
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Appendix  C.  Discriminant  Analysis  Results  (42  cases) 


C-l.  Classif icat ion  based  on  C413  w/o  PTIME 


Funct ion 


Eigenvalue 

.165 


After  Funct ion  W i 1 k ' s A  Sig 
0  .858  .051 

1  - . - 


Standardized  Discrimination  Coefficients 


Fund 


TIME 


1.000 
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Classification  Results 

Predicted 

Group 

Membership 

* . ' .  * 

LV-V- 

•V\> 

i  ’  « 

Actual  Group 

#  of  cases 

1 

2 

3 

■  v-v- 

Group  1 

18 

4 

(22.2) 

5 

(27.8) 

9 

(50.0) 

r  '-Ta 
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•.•••->. 

Group  2 

20 

1 

(  5.0) 

11 

(55.0) 

8 

(40.0) 

rV.V, 

Group  3 

4 

0 

(  0.0) 

1 

(25.0) 

3 

(75.0) 

%  of 

grouped  cases 

correctly  classified  - 

42.86 

C-2.  Class  if icat ion  based  on  CP2 


After  Function 

Wilk'sX  Sig 

Funct ion 

E  igenvalue 

0 

.700  .035 

1 

.315 

1 

.921  .208 

2 

.086 

Standardized 

Discr iminat ion 

Coeff icients 

Func  1 

Func2 

TIME 

-.64768 

.41345 

EAVG 

.14314 

.89772 

DT 

.70203 

.50386 

Classification  Results 

Actual  Group 

#  of  cases 

Predicted  Group 
1  2 

Membership 

3 

Group  1 

17 

10  5 

(58.8)  (29.4) 

2 

(11.8) 

Group  2 

13 

4  8 

(30.8)  (61.5) 

1 

(  7.7) 

Group  3 

12 

4  4 

(33.3)  (33.3) 

4 

(33.3) 

%  of  grouped  cases  correctly  classified  -  52.38 


C~3.  Class  if icat ion  based  on  CS1  w/o  PTIME 


After  Funct ion  Wilk ' s>  Sig 


Funct ion 

Eigenvalue 

0 

.798  .069 

1 

.  149 

i 
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2 
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Standardized 
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1.00589 

-  .32728 
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3 
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2 
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Group  3 

16 

4  7 

(25.0)  (43.8) 

5 

(31.3) 

%  of  grouped  cases  correctly  classified  -  47.62 


C-4.  Class  if icat ion  based  on  CS1  w/PTIME 


After  Funct ion  Wilk 1 s  A  Sig 


Funct ion 

Eigenvalue 

0 

.749  .089 

1 

.164 
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.872  .074 
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Discrimination 
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Func2 
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PTIME 
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Classification  Results 

Actual  Group 

#  of  cases 

Predicted  Group 

1  2 

Membership 

3 

Group  1 

12 

4  5 

(33.3)  (41.7) 

3 

(25.0) 

Group  2 

14 

1  9 

(  7.1)  (64.3) 

4 

(28.6) 

Group  3 

16 

4  7 

(25.0)  (43.8) 

5 
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%  of  grouped  cases  correctly  classified  -  42.86 


Appendix  D.  Discriminant  Analysis  Results 
(T=l ,  30-38  cases) 

D-l.  Classif icat ion  based  on  C413  (32  cases) 


After  Funct ion  Wilk ' s  X  Sig 


Funct ion 

Eigenvalue 

0 

.679  .026 

1 
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Standardized 

Discrimination 
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%  of  grouped  cases  correctly  classified  -  53.13 


D-2.  Classif icat ion  based  on  CPRMS  (33  cases) 


After  Funct ion  Wilk 1 s X  Sig 
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4 
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%  of  grouped  cases  correctly  classified  -  48.48 
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